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Abstract. A review of the muon spin rotation and relaxation (/xSR) studies on 
magnetic materials published from July 1993 is presented. It covers the investigation 
of magnetic phase diagrams, of spin dynamics and the analysis of the magnetic 
properties of superconductors. We have chosen to focus on selected experimental 
^ ■ works in these different topics. In addition, a list of published works is provided. 
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1. Introduction 

The breadth, format and style of this article is intended to provide an accessible 
and stimulating review of recent investigations of the physical properties of magnetic 
materials by the fiSK experimental technique. /iSR is an acronym for Muon Spin 
Rotation, Relaxation or Resonance. We shall only deal with the first two of these 
techniques, i.e. the two most commonly used in studies of magnetic materials. 
Measurements can be carried out with positive and negative muons. Since almost 
all the investigations were performed using the positive muon, we shall only discuss 
these. By no means have we attempted to write a comprehensive review. At the risk 
of being invidious, we discuss in three sections selected recent works which display 
the possibilities of the technique. To achieve a balanced picture, we list the works 
published from July 1993. The material published before July 1993, including that 
presented at the /iSR conference held in early 1993 on the island of Maui, Hawai, has 
already been nicely reviewed (Schenck and Gygax 1995). 

The organisation of this article is as follows. In section ^ we introduce the basic 
concepts of /xSR. In section || we discuss magnetic phase diagram studies. Section || 
presents two examples of the investigation of spin dynamics in magnets. The next 
section (section ||) deals with the very successful studies of the magnetic properties of 
superconductors. In section || we summarize the present status of /iSR and mention the 
scheduled technical developments at the /iSR facilities. In the last section (section [?]) 
we give a list of published works. This review is completed by four appendices, the 
material of which is partly original. 

2. /iSR: Muon Spin Rotation, Relaxation 

Since the /xSR technique has been described in many reports (Seeger 1978, Chappert 
and Grynszpan 1984, Schenck 1985, Chappert and Yaouanc 1986, Cox 1987, Karlsson 
1995, Schenck and Gygax 1995, Schatz and Weidinger 1995), we will only sketch it 
briefly. More information is provided in the appendices. 

Currently, /iSR experiments can be performed at three facilities located at i) 
TRIUMF (4004 Wesbrook Mall, Vancouver BC, Canada V6T 2A3), ii) the Paul 
Scherrer Institut (PSI /xSR Facility, CH-5232 Villigen PSI, Switzerland) and iii) the 
Rutherford Appleton Laboratory (ISIS Facility, Chilton, Didcot, Oxon 0X11 OQX, 
United Kingdom). In addition, measurements are carried out at the Meson Science 
Laboratory of the Faculty of Science of the University of Tokyo (Bunkyo-ku, Tokyo 
113, Japan) and the Phasotron of the Joint Institute for Nuclear Research (Dubna, 
Head Post Office, P.O. Box 79, Moscow, Russia). The muon beams at the ISIS 
Facility and Meson Science Laboratory being pulsed beams are well adapted to study 
weak magnetic signals: this capability results from the virtual absence of background 
related to contamination of the beam with particles other than muons. However, 
their relatively low time-resolution means they are unsuitable for investigating systems 
exhibiting fast relaxation processes or appreciable local fields (larger than ~ 50 mT) 
at the muon site. The latter limitations do not apply to the beams at the other 
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institutions since they are quasi-continuous. In fact, the muon beams at the different 
institutions are complementary. 

The /zSR technique uses the positive muon as a probe. The muon may form a 
bound state with an electron, called muonium, an exotic isotope of hydrogen. However 
muonium has never been observed in metals and the probability of its formation in 
oxides is small. We shall therefore consider only free muons. 

In the /iSR technique polarized muons are implanted into a sample where their 
polarization evolves in the local magnetic field until they decay (the muon lifetime 
is 2.2 us). Because of its positive charge, the muon localizes at an interstitial site. 
Due to the absence of quadrupolar electric moment (spin 1/2) the muon does not 
couple to electric field gradients. The decay positron is emitted preferentially along 
the muon spin direction; by collecting several million positrons, one can reconstruct 
the time dependence of the muon spin-depolarization function which, in turn, reflects 
the spatial and temporal distribution of magnetic fields at the muon site. Because of 
their large kinetic energy (« 30 MeV), the positrons easily go through the sample. 
Fortunately they are weakly absorbed by cryostat or furnace walls so that complex 
sample-environment equipment may be used. 

Two types of experimental geometries are generally used, see figure @. In the 
longitudinal geometry an external magnetic field B cxt is applied along the initial 
muon beam polarization direction S M and positron detectors are set parallel and 
antiparallel to S M . We then refer to the forward and backward direction, respectively. 
In the transverse geometry B ext is perpendicular to S M and the positrons are 
detected perpendicular to B cxt . The zero- field measurements are performed with 
the longitudinal geometry. 

The /xSR technique is a local-probe hyperfine-method as are nuclear magnetic 
resonance (NMR), perturbed angular correlations (PAC) and Mossbauer spectroscopy. 
Therefore the basic physical concepts developped for the latter techniques can be 
transferred to /iSR. One of its unique characteristics is its ability to detect very small 
magnetic moments. This is understood as follows. 

We first recall that a magnetic moment precesses around a local magnetic field 
Bi oc at a pulsation frequency to proportional to B\ oc : this is the Larmor precession. 
For the muon u)^ — 2ixv ^ = 7p-Bi oc where is the precession frequency and 7^ is 
the muon gyromagnetic ratio (7^ = 851.6 Mrad s _1 T _1 ). Therefore, in a constant 
field, the moment rotates by an angle 7^-Bioci in the elapsed time t. Since an angle 
of 0.5 radian is measurable (cos 0.5 ~ 0.88) and a /iSR measurement can be routinely 
carried out up to t = 15 /its, a local field as small as 0.04 mT can be detected. It 
might be produced by a nuclear magnetic moment which is ~ 10 3 times smaller than 
an electronic magnetic moment. 

While the /iSR technique can detect very small magnetic fields, very high fields 
can not be measured: with a conventional spectrometer of a continuous muon source 
the maximum field is ~ 3 T. Recent measurements at 6.5 T have been performed by 
Riseman et al 1995. 

The muon spin depolarization function is written as aP a (t). It is also called 
the asymmetry, a is the initial asymmetry (at t = 0) and P a (t) the normalized 
depolarization function which will be referred to as the depolarization function. The 
Cartesian label a denotes the direction along which the muon polarization is measured, 
i.e. it can be X, Y or Z. The value of a depends on the experimental geometry; 
typically a — 0.25. We use an orthonormal laboratory reference frame. B cxt is taken 
along the Z axis and parallel to Z or X in the longitudinal and transverse set- 
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Figure 1. Principle of the two types of experimental geometries: the transverse 
and longitudinal set-ups arc cartooned in the upper and lower panel, respectively. In 
order to produce drawings which are easy to understand, the muon beam momentum 
and the initial muon beam polarization have been taken parallel. In reality, these 
two vectors are antiparallcl. The arrows originating from the sample sketch the muon 
spin direction: the dashed one at t = and the solid one at t. The dashed and solid 
line cardioids drawn for the transverse field setup panel represent the probability 
of positron emission in a given direction relative to the muon spin position at the 
instant of the decay. 
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up, respectively. Whereas in a longitudinal set-up, only Pz(t) is of interest, in the 
transverse set-up, Px{t) and Py(t) can be measured. Since these two functions always 
contain the same information, we only consider Px(t). 

P a (t) monitors the properties of the magnetic field at the muon site. If all the 
muon spins precess in the same static magnetic field, oriented at an angle 6 from S M , 
the Larmor equation yields 



This is the result on which the entire fiSH technique is based. 

We consider a magnet in zero external magnetic field. Then Pz(t) is of interest. 
Bi oc is usually not zero below the magnetic phase transition. For a magnet in 
polycrystalline form, the spatial average of (1) has to be performed. If the sample 
is not textured, we obtain 



The oscillating component reflects the magnetic order in the sample. This single 
component exists even in a polycrystalline sample because we suppose there is only 
one type of muon localization sites and for these sites Bi oc is the same. If the magnet 
is disordered, i.e., if its correlation length is small, Bi oc can take a large number 
of values. Then the oscillation can be strongly damped and even disappear. If the 
muon spins precess too quickly relative to the time resolution of the spectrometer, 
the oscillation will not be observed and the resulting muon spin-polarization will be 
averaged out to 0. Then, at a magnetic phase transition, if no wiggles are observed in 
the muon signal, one expects a drop in the effective initial asymmetry from a in the 
paramagnetic state to a/3 in the ordered state. In figure |^ we present examples for 
two samples. 

Whereas for GdNis we do observe the expected behaviour, for UPt2Si2 the change 
at T/v is smoother than expected and the ratio of the effective initial asymmetry in 
the paramagnetic and ordered state is smaller. We then conclude that the distribution 
of Tn values is relatively large in the UPt2Si2 sample. In addition, this sample has a 
strong texture. The distribution is the signature of crystalline disorder: probably some 
Pt and Si atoms interchange their atomic positions. Complementary measurements 
at ISIS on a single crystal sample of UPt2Si2 with perpendicular to Bi oc (and 
therefore to the c axis) show that the effective initial asymmetry decreases to zero 
around T/v in a temperature interval as large as ~ 5 K (Gubbens et al 1996). This 
result supports the conclusion deduced from the polycrystalline sample measurements. 
Note that the measurement of the effective initial asymmetry offers a stringent test of 
the sample quality. 

Figure |2| shows that the maximum initial asymmetry is smaller than expected. 
The difference is accounted for by the contribution from the background, i.e. muons 
which have not been stopped in the sample, but for instance in the sample holder. We 
write aPz(t) = a s P|(t) + abg-P^ s (i) where the first and the second term describes 
the contribution from the sample and the background, respectively. a s is plotted in 
figure |[ Usually the sample holder is a silver plate, P^if) is then taken as time 
independent (P^ s (t) — 1) because silver has no electronic moment and very small 
nuclear moments. On the other hand, if one analyses an extremely weakly damped 
/iSR signal, the damping due to the silver nuclear moments should be taken into 
account. 



P a (t) = cos 2 9 + sin 2 0cos(u; M t). 



(1) 



Pz(t) = ~ + ~cos(uV). 



(2) 



() 



Dalmas de Reotier and Yaouanc: fig. 2a 



A- 



GdNi 5 

zero-field 



(a) 



20 40 60 ; 
Temperature (K) 



100 



a> 
E 
E 



LU 
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Figure 2. Temperature dependence of the effective initial asymmetry after 
backgroung substraction measured in two polycrystalline samples at the ISIS Facility 
(Dalmas de Reotier et al 1990 and Gubbens et al 1996 ). Whereas this asymmetry 
drops at Tq in a very small temperature range (~ 0.1 K) in the ferromagnet GdNis 

(a) , it changes gradually near Tjv (within ~ 5 K) in the antiferromagnet UPt2Si2 

(b) . The solids lines are guide to the eyes. 



The depolarization of the muon spins is caused either by a static distribution 
of local fields (dephasing of the muon spins) or by fluctuating fields (relaxation of 
the muon spins). The general framework needed to describe the combination of the 
two mechanisms will be first given. In section 2.1 and section 2.2, respectively, the 
resulting Pz(t) and Px(t) will be presented. The relaxation induced by a purely 
dynamical process described in section 2.2. 

When measuring P a (t), an ensemble average is made. Usually the muons do 
not all experience the same local field, i.e. there is a field distribution at the muon 
site, .D(Bi oc ). The time evolution of the muon magnetic moment experiencing a field 
Bi oc can be computed by solving the Larmor equation in classical mechanics (or the 
Heisenberg equation in quantum mechanics). We define S^ a (t, Bi oc ), the projection 
along the a axis of the unit vector parallel to the muon moment. The depolarization 
function is simply written 



Pa(t) — J Sn itx (t, Bio C )D(Bi oc ) cEBioc 



(3) 



If the field distribution is sharp, _D(Bi oc ) is a Dirac delta function and we recover 
(1). Different origins for £>(Bi oc ) exist depending on the physical case. For a 
magnet the distribution is due to the nuclear and electronic magnetic moments. For 
a superconductor in the mixed phase, the nuclear magnetic moments and the flux 
quanta induce the distribution. First, for simplicity, we suppose that this distribution 
is static, isotropic and Gaussian centered at 0. Therefore the distribution is written 
as D(B loc ) = D(B* c ,BZ c ,B? oc ) where 
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A 2 '/if,, is the variance of the components of the field distribution. A can be expressed 
in terms of the interaction between the muon spin and the spins of the compound and 
the characteristics of the latter spins (see Schenck 1985 for details). In contrast to 
NMR where the field at the probe of dipolar origin is negligible (except in the case of 
NMR on protrons), the dipolar interaction dominates in the /iSR case. This can have 



important concequences as shown for example in section 4.1. 

Below, we shall present the depolarization function for the two possible 
experimental geometries, and discuss the meaning of the magnetic field at the muon 
site. 

2.1. The transverse depolarization function 

Up to the end of the eighties most of the experiments were performed with the 
transverse geometry (or in zero external field for the case of magnetically ordered 
materials), mainly because the spectra are then less sensitive to the quality of the 
muon beam, e.g., contamination with other particles, than when recorded with 
the longitudinal geometry. Today the transverse geometry is mainly used for the 
investigation of the microscopic field distributions arising from either the vortex state 
in superconductors or magnetic domains in metals. In addition, a very active field 
of investigation is the measurement of the local susceptibility at the muon site. This 
type of experiment yields information on the symmetry at the muon site and on the 
hyperfine coupling-constant. 

The field distribution leads to a dephasing of the muon spins responsible for the 
depolarization. If _B oxt is sufficiently large relative to i.e. if B cxt > 5(A/j ll ) 

(Dalmas de Reotier and Yaouanc 1992), Px(t) probes only the field distribution along 
B cxt , i.e. along the Z axis. Then, using (3) and (4), the following simple result is 
found: 

/ A 2 t 2 \ 

Px(*) = expl — J cos(w M t). (5) 

If the argument of the Gaussian term is small, i.e. if t A is small, the envelope of Px if) 
is well approximated by the parabolic form (1 — A 2 t 2 /2). 

Usually the field distribution is not static. A useful first approximation is to 
account for the dynamics with a single fluctuation rate, and to use the stochastic 
theory of dynamical processes. Kehr et al 1978 have used the strong collision model 
which supposes that B\ oc takes a given value for a time 1/v followed by a new value 
not related to the previous one, i.e. they consider a Markov process. They have shown 
that the following analytical formula is a fair approximation: 

Px(t) = cxp j-^- [exp(-i/t) - 1 + i4] J cos(w„t). (6) 

This formula, called the Abragam formula (in fact it was first derived by Anderson 
1954 as a model for the NMR line shape), interpolates between the static case, (5), 
and the fast fluctuation limit, the so-called motional-narrowing limit, for which the 
envelope is an exponential function: 

P x (t) = exp (-Ax*) cos(ovi), A x = A 2 jv. (7) 

Although the Gaussian field distribution model with the dynamics described as a 
Markov process provides only a rough picture of the physics involved, it is useful since 



8 



it clearly shows the physical origin of the measured damped oscillation. The limits of 
this model for the depolarization induced by nuclear magnetic moments is discussed 
by Dalmas de Reotier et al 1992, Yaouanc and Dalmas de Reotier 1994 and Cameron 
and Sholl 1994. In a magnet, the functional form of D(Bi oc ) reflects its magnetic 



structure. This is discussed in Appendix A. In a type II superconductor, _D(Bi oc ) is 



usually not Gaussian. In fact its shape is characteristic of the vortex state as shown 
in section [|. Additional information is provided in Appendix E . 



2.2. The longitudinal depolarization function 

Nowadays, the longitudinal geometry is popular for the study of magnets since it is 
well suited for the characterization of a magnetic phase transition. The fact that the 
measurements can be performed in purely zero-field is a definitive advantage of the 
/iSR method. 

In zero-field, an analytical formula is found for Pz(t) if the distribution is isotropic 
and Gaussian. Using (3) and (4) one derives 

Pz{t) = \ + \(l-AH*)eM-\^t 2 )- (8) 

This is the well-known Kubo-Toyabe function which is plotted in figure ||. Pz(t) 
exhibits a dip at t — y/3/A. If tA is sufficiently large, it saturates to 1/3. If tA is 
small, Pz(t) is well approximated by the parabolic form Pz(t) — 1 — A 2 t 2 . Relative to 
the transverse case, the initial depolarization is stronger in zero-field by a factor of two 
since for the latter geometry the two components of Bi oc perpendicular in the Z axis 
participate to the depolarization. Equation (8) is strictly valid for an isotropic field 
distribution at the muon site. Szeto 1987, Ealmas de Reotier 1990 and Solt 1995 have 
computed Pzif) for different types of anisotropic field distributions. As expected, the 
value of Pz(t) for tA large reflects the type of anisotropy. Equation (8) is derived 
with the hypothesis that the average field at the muon site is zero. An extension for 
a finite average field inside a polycrystalline sample has been given by Kornilov and 
Pomjakushin 1991. 

As for the transverse geometry, the dynamics can be accounted for approximately 
by the strong-collision model (Hayano et al 1979). Examples are shown in figure [| 
Interestingly, at large t, the functional form of Pz(t) depends strongly on v when the 
distribution is quasi-static, and the Kubo-Toyabe function can still be used but with 
the 1/3 factor multiplied by exp(— 2i4/3). Therefore, slow dynamical processes can 
be characterized by zero- field measurements. If v j A is sufficiently large, Pz(t) is well 
approximated by the analytical formula (Dalmas de Reotier and Yaouanc 1992) 

f 2A 2 1 

P z (t) = cxp <^ -— [exp(-ut) -l + vt]\. (9) 

This approximation is compared in figure ^| to the exact numerical solution for 
v/A = 3. This analytical formula is the envelope of the transverse depolarization 
function (see (6)), except for the substitution A 2 — > 2A 2 . The factor 2 reflects 
again the fact that the depolarization is induced by the two components of Bi oc 
perpendicular to the Z axis. Further analysis of Pz (t) in a low field for time-dependent 
random magnetic fields is given by Shibata and Shimoo 1995. In the motional 
narrowing limit, i.e. when the dynamics is fast, v -C 1/t, Pz(t) is an exponential 
function 

Pz = exp(-Azt), X z = 2A 2 /u. (10) 
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Figure 3. (a) Fluctuation-rate dependence of the zero-field depolarization function 
derived for a isotropic Gaussian field distribution and a Markov dynamical process, 
(b) Magnetic field dependence of the longitudinal depolarization function derived for 
a static isotropic Gaussian field distribution. The values of the fluctuation rate and 
magnetic field are respectively given in units of A and The dashed line in 

the (a) panel is the prediction of the modified Abragam formula ((9)) for v/A = 3. 
The Kubo-Toyabe function (Eq. 8) corresponds to the curves marked with in both 
panels. 



In order to fully characterize a quasi-static process, it is important to perform 
longitudinal field measurements. Figure [3] presents the field dependence of Pz(t) 
for a static field distribution: Pz{t) is strongly field dependent. Obviously, the strong 
collision model can be applied to this depolarization functions to account for dynamics. 
Recently Uemura et al 1994 have pointed out that the following scaling property holds: 

P Z (A, Bext, v, ft) = P z (fA, /B cxt , fv, t). (11) 

Physically, / can be the fraction (/ < 1) of the time during which the muon experiences 
a non-zero magnetic field, i.e. in the time fraction (1 — /), this field is zero. Spin liquid 
systems could be an example of such a model (Uemura et al 1994, Bonville et al 1997). 
Note that fitting a spectrum with a regular dynamical Kubo-Toyabe function with a 
longitudinal field then yields reduced effective /A, fB cyLt and fv values. 

The field distribution at the muon site in disordered systems is Gaussian only 
when the spin concentration is large. Walstedt and Walker 1974 have shown that 
the distribution is Lorentzian for a disordered dilute spin system. Kubo 1981 has 
computed Pz (t) for such a distribution. In the motional narrowing limit Pz (t) is then 
a square root stretched exponential function. Pz{t) for disordered systems has been 
discussed recently by Berzin et al 1993, Borgs et al 1995 and Crook and Cywinski 
1997. 

Even if the Gaussian approximation is fair, as it is usually the case for the 
depolarization due to nuclear spins, the Kubo-Toyabe formula and its extension 
to dynamical processes does not provide a sufficiently good model to analyse high 
statistics spectra. The interaction between the muon spin and the lattice spins has to 
be taken into account, at least in nuclear systems (Celio 1986, Dalmas de Reotier and 
Yaouanc 1992, Keren 1994, Yaouanc and Dalmas de Reotier 1995). 

For a purely fast dynamical process, the relaxation of the muon spin leads to 
an exponential depolarization function characterized by a relaxation rate \z\ see 



10 



Appendix B. This process is very similar to the NMR spin-lattice relaxation, and 
consists of an exchange of energy between the lattice spins and the two muon 
Zeeman levels. In that case, the name of relaxation function for Pz(t) seems more 
appropriate. We point out that Pz(t) is an exponential function only if no spatial 
average is necessary. Therefore, in general, Pz{t) is not an exponential function for 



a polycrystalline sample (see Bonville et al 1996). In Appendix C we discuss the 



meaning of Xz in terms of the spin correlation-functions of the magnet. 
2.3. The magnetic field at the muon site 

A detailed discussion of the field at the muon site, Bi oc , is given by Schenck and Gygax 
1995. 

In general Bi oc , is the sum of seven terms: 

Bloc = B con + Btrans + B dip + Bl + B d em + Bdia + B cxt . (12) 

B C on is the contact hyperfine-field resulting from the spin density at the muon site 
which is induced by the polarization of the conduction electrons. Therefore it only 
exists in metals. Btrans is the transferrred hyperfine field. In metals this field is 
due to the indirect Rudermann-Kittel-Kasuya-Yosida (RKKY) interaction. The next 
three terms in (12) reflect the muon spin interaction with the localized lattice spins 
through the dipolar interaction. This interaction gives rise to a dipolar field which 
is expressed as a lattice sum over the sample. This sum is split into two parts by 
separating the volume of the sample into a sphere around the muon (the Lorentz 
sph ere) and the rest. B dip is given by a lattice sum restricted to the Lorentz sphere; 



see Appendix A. Summing over the rest yields the Lorentz and the demagnetization 
fields, B L (= (jio/3)Msat) and B dcm (= -/i NM bu lk), respectively. M sat and M bu ik 
are respectively the saturation and bulk magnetizations, /Lto the permeability of free 
space and N the demagnetization factor tensor which depends only on the sample 
shape. The trace of N is 1. For a sphere Bdcm = — MoM DU ik/3. For an infinite plane, 
in practice for a disk with extremely small thickness relative to radius, B d cm = 
if B ox t is applied perpendicular to the disk axis and Bdcm — — /^o-M-buik if B C xt is 
parallel to the disk axis. Bdia in (12) is the diamagnetic field which is important only 
for superconductors. In usual magnets, it can be neglected. 

In a magnetized ferromagnet all the terms of (12), except Bdia, are to be taken 
into account. In a ferromagnet which is not macroscopically magnetized, Bd cm = 
because Mbulk = 0- For an antiferromagnet in zero field, Bl and Bdcm are necessarily 
zero. In a metal, B con and B tra ns are normally independent of the crystal direction. 
B dip depends strongly on the symmetry at the muon site; for examples, see Seeger 
1978. 

For a paramagnet or a superconductor in an external field, a shift of the muon 
frequency is usually observed, i.e. Bi oc is different from B ox t. The induced field 
(Bioc — B xt) is not necessarily parallel to B GX t- However, in general, |Bi oc — B ex t| 
•C | B cx t | . Therefore it is useful to characterize the observed shift by the projection of 
(Bioc — Bcxt) onto Bcxt: 



B ex t • (Bi oc — B C xt) 

3 cxt 



^cxp _ -"cxi v-"loc "cxt; ^g-j 



Since K cxp is not intrinsic to the compound (it depends on the sample shape), one 
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defines a new ratio, K, traditionally called the Knight shift 

B cx t • ^B con + Btrans + Bj ip ^ 

-°ext 

The two ratios can be related using (12) 



> m ^- (14) 



K = K cxp - Q - n) XP ~ (1 - n) Xd , (15) 

where \p is the paramagnetic susceptibility, Xd the diamagnetic susceptibility 
produced by persistent currents in the mixed state of a superconductor and n the 
appropriate component of N (Heffner et al 1989). Obviously Xd = if the compound 
is not superconducting at the temperature of the measurement. At first sight, one may 
be surprised that X p an d Xd are not multiplied by the same factor. This is understood 
if one remembers that Bl responsible for the 1/3 factor originates from the dipole 
moments in the Lorentz sphere. This concept is unsuitable for the superconducting 
electrons. The factor 1 in (1 — n) originates from Bdi a , i.e. it does not concern the 
paramagnetic properties of the compound. 

The concept of Knight shift is useful because one expects this quantity to be 
field independent. In particular, in superconductors, its thermal behaviour when 
crossing the superconducting temperature yields direct information on the parity of 
the superconducting order parameter. However, recent measurements on UBei3 and 
UBe 12 .gi-Bo.o9 indicate that the Knight shift concept is not valid for these compounds 
since \B\ oc {T) — B\ oc (T c )} is field independent (Luke et al 1991, Heffner et al 1997), 
i.e. B\ oc {T) does not scale with B ext . 

Since K depends linearly on the components of the susceptibitity tensor, for a 
given B cxt direction, it is presented as a function of the corresponding susceptibility, 
the temperature being an implicit parameter. From such plots (the Clogston-Jaccarino 
plots) the value of the components of the dipolar tensor and of the hyperfine coupling 
constant can be determined (Schenck and Gygax 1995). 

The study of the angular dependence of the local field at the muon site measured 
in transverse field on a single crystal in the paramagnetic phase gives the opportunity 
of determining the muon site. In figure Qjsuch an angular dependence is shown for 
the cubic system CeBg (space group PrnZm). The presence of two different fields 
(as seen in the Fourier transform of the mesured signal) indicates the existence of 
two magnetically inequivalent sites for the muon. From the analysis of this angular 
dependence and the value of the component of the dipolar field deduced from the 
Clogston-Jaccarino plot, Amato et al 1997 infer that the muon occupies the d site 
(0,0,1/2). In recent years this type of procedure has been used to determine the 
muon site in many intermetallic compounds. The results show that it is not possible 
to predict reliably the muon site localization. For a given crystal structure, this 
site depends on the chemical formula of the compound. It may even depend on the 
temperature. Because of its electric charge, the muon may distort locally the crystal 
lattice and induce an electric field gradient on its neighbor atoms. 



3. Magnetic phase diagrams 



The /xSR technique has been popular in the condensed matter community for its 
success in detecting magnetic phase transitions in compounds with small magnetic 
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Figure 4. Angular dependence of the two local fields observed in CeB6. B ox t 
was rotated around the [1,1,0] axis. The two signals have an amplitude ratio 2:1 
corresponding on one hand to the (1/2, 0, 0) and (0, 1/2, 0) sites (open symbols) and 
on the other hand to the (0, 0, 1/2) site (solid symbols). The solid lines are fit to the 
data (adapted from Amato et al 1997). 



moments. In this section we examine four such examples. In section 3.1 and section 3.2 
we discuss recent results obtained for the organic and low dimensional magnets, 
respectively. Whereas for the organic compounds oscillating /iSR signals are detected, 
pointing to the existence of relatively well ordered magnetic strutures, in the 3-leg 
spin la dder s one only observes wide field distributions, a fingerprint of disorder. In 
section 3J3 we show that the /xSR technique can be used to detect magnetic domains 

The appearance of the domains is dir ectly 
Finally we present in section 3.4 new 



in a diamagnetic metal such as beryllium, 
connected to the de Haas- van Alphen effect, 
results obtained for some heavy fermion metals in which two / electron components 
exist. 



3.1, Magnetic ordering in organic compounds 

The search for purely organic molecular ferromagnets which contain only light elements 
(carbon, nitrogen, hydrogen and oxygen) is a subject of strong current interest. The 
first such material to be found, the f3 crystal phase of para-nitrophenyl nitronyl 
nitroxide (p-NPNN, Ci3Hi 6 N 3 04), was reported to have a Curie temperature Tc 
~ 0.6 K (Tamura et al 1991). The unpaired spin is associated with the nitronyl 
nitroxide group (N-0 group). The role of the rest of the molecule is to ensure the 
appropriate overlap of the correct orbitals on neighbouring molecules to produce 3D 
ferromagnetism. A whole series of materials which incorporate this N-0 group has 
been synthesized. 

The first direct observation of spontaneous magnetic order in p-NPNN was done 
using the fiSR technique (Le et al 1993b). It has been subsequently confirmed by zero- 
field neutron diffraction (Zheludev et al 1994). In figure ||a we present three zero-field 
fiSK spectra. The mere observation of an oscillating signal at low temperature, i.e. of 
a spontaneous internal magnetic field B\ oc , is a clear signature of the existence of static 
magnetic correlations. Since a long-lived oscillation is detected rather than an increase 
in damping, p-NPNN orders with a well defined magnetic structure. Measurements 
on crystals indicate that Bi oc , is nearly parallel to the b crystal axis which is the easy 
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axis (Le et al 1993b). The shape of the spectra changes drastically between 650 mK 
and 700 mK as seen in figure the magnetic phase transition occurs between these 
temperatures. In figure [^b we present the temperature dependence of B\ oc . The solid 
line is a fit with B\ oc {T) oc [1 — (T/Tc) Q ] • This compact formula allows us to discuss 
the spin wave (T < T c ) and the critical regimes: for T < T c , [Bi oc (0) ~ Bi oc (T)} oc T a 
and Bioc(T) oc {T-T c f near Tc- B\ oc (T) is well described in the whole ferromagnetic 
state by the compact formula with a = 1.7 (4) and (3 = 0.36 (5) (Le et al 1993b and 
Blundell et al 1995). These results are consistent with that of a 3D Heisenberg magnet. 
The weak magnetic anisotropy is accounted for by the dipolar interaction between the 
unpaired electron spins (Le et al 1993b). The small B\ oc {T = OK) value found in 
the organic magnets is a strong indication that the magnetic moment carried by the 
unpaired electrons is small since it would be surprising that the dipolar field at the 
muon site accidentally cancels almost perfectly in all these compounds. An analysis 
of the behaviour of the oscillations in an applied longitudinal field, in terms of the 
magnetization process and demagnetising field, is given by Blundell et al 1995. 
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Figure 5. (a): some zero-field spectra recorded for p-NPNN (adapted from Blundell 
et al 1995). The solid lines are fits, (b): temperature dependence of the local 
magnetic field at the muon, B loc , in p-NPNN, p-PYNN and 3-QNNN (adapted from 
Blundell et al 1994 and 1995). The solid line is a fit with B loc (T) oc [1 - (T/T c ) a ]' 3 . 
The dashed lines are fits with the molecular field model with spin S = 1/2. 



The transition temperature and the nature of the magnetic ground state of the 
molecular crystals based on the nitronyl nitroxide radical depend strongly on their 
crystal structure. This is nicely seen in figure §3: T c ~ 90 mK for p-PYNN , T c ~ 
230 mK for 3-QNNN and T c w 670 mK for p-NPNN. Other examples and a detailed 
discussion are given by Blundell et al 1997b. 

It is well known that neutron scattering can provide extensive information on 
magnets, in particular on the magnetic structure and on the excitations. But in 
compounds with small magnetic moments such as the organic magnets, the ^SR 
technique is a sensitive and useful probe: it easily yields the value of the critical 
temperature and information on the thermal behaviour of the order parameter. In 
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addition, it should be possible to investigate the excitations at a very small energy 
transfer (see section 4.1 for an example for a conventional ferromagnet). 



3.2. Magnetic ordering in spin chains and ladders 

Hoping to gain insight into the mechanism of superconductivity, some experimentalists 
and theoreticians have recently turned their attention to systems as simple as one 
dimensional chains and ladders of copper and oxygen atoms with the aim to eventually 
apply their findings to the high T c superconductors where Cu-0 planes play an 
important role. The //SR technique has been used to determine the magnetic 
properties of the ground state of some of these systems. 

The ground state properties of a linear chain of antiferromagnetically coupled 
spins have been intensively studied because of the pronounced quantum effects. Both 
integer and half integer spin chain systems have a singlet ground state (Mermin and 
Wagner 1966, Haldane 1983a and 1983b). Because the spin excitations are gapless 
at momentum q — and 7r for half integer spin chains (des Cloizeaux and Pearson 
1962), a magnetic ordering can be expected when interchain interactions are taken into 
account. However, for integer spin chains which have a so-called Haldane energy gap 
(Haldane 1983a and 1983b), no ordering should occur. The key parameters for half 
integer spin chains are the magnitude of the ratio Tjv/J {J is the intra-chain coupling 
constant) and the ordered magnetic moment at T — K, M(T = 0). J is estimated 
from magnetic susceptibility and infrared light absorption measurements. The Neel 
temperature is easily determined by zero-field /iSR in such systems characterized by 
small magnetic moments. Neutron diffraction is the most direct method to measure 
M (T = 0). Although the ^SR technique cannot yield a precise M(T = 0) value for a 
given compound, it is well suited to measure accurately the relative size of moments 
for iso-structural materials. Therefore a combined neutron diffraction and //SR study 
is expected to yield reliable results. Kojima et al 1997 have done such a study 
for the quasi one-dimensional antiferromagnets Sr 2 Cu03 and Ca2CuC>3. Plotting 
M(T = 0) versus Tn / J and comparing with different models, they find that the 
chain moan field approach best explains the experimental results. This approach takes 
more quantum effects into account than the two different spin-wave approximations 
available. Probably, this is the reason why it provides a better description since the 
moment reduction is dominated by quantum spin fluctuations. 

As mentioned, the spin ladder systems have been studied because it is thought that 
an understanding of their magnetic properties is a prerequisite for a proper description 
of the magnetic properties of infinite CuO layer systems. A 3- leg ladder structure is 
displayed in figure |g. The oxides Sr„_ 1 Cu„ + i02„ are realizations of such ladders. 
Indeed, one observes that the geometry of the ladder structure and of the CuO square 
lattice layer are related. A nice review of the physics of these systems is given by Goss 
Levi 1996. Note that the physics of the spin ladder systems and of the Haldane spin 
chains are closely related (Strong 1997). 

A key prediction is that only ladders with even numbers of legs have a singlet 
ground state separated from the triplet state by a large spin gap (Rice et al 
1993). However, the odd- leg systems are expected to reach a magnetically ordered 
ground state in the presence of interladder interactions. Kojima et al 1995a have 
performed zero-field and longitudinal field /iSR measurements to test these theoretical 
predictions. 

In figure u\ spectra recorded on the 3-leg system are presented. The strong 
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Dalmas de Reotier and Yaouanc: fig. 6 
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Figure 6. The 3-leg ladder structure (from Kojima et al 1995a). Oxygen ions locate 
at each corner of the squares. The ferromagnetic interladder interaction J' is much 
smaller than the antiferromagnetic intraladdcr interaction J. 



depolarization of the zero-field spectra at low temperature shows that the ground 
state is magnetic. Since no wigg les a re detected (in contrast to the observations for 
the organic magnets; see section 3.1), the disorder in the compound is important or 
the number of muon localization sites with different local fields is large. Comparing 
the spectra recorded at 50 K and 60 K, we infer that a 3D magnetic phase transition 
occurs between these temperatures. The longitudinal field measurements confirm the 
interpretation of the zero-field spectra, i.e. the ground state of the 3-leg system with 
interladder interactions is a conventional static ordered state rather than a spin liquid 
system since / — 1 (see (11)). 
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Figure 7. Some spectra recorded on Sr4Cu60io which has a 3-leg spin ladder 
structure. The solid lines are fits (adapted from Kojima et al 1995a). 



The magnetic behavior of the 2-leg and 3-leg ladder systems differ remarkably 
as seen in figure pi The depolarization functions are described with a square-root 



exponential function, appropriate for dilute fluctuating moments (see section 2.2). 
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Therefore no static magnetic ordering is detected. The depolarization originates from 
dilute unpaired spins which may be associated with defects in the sample. 

In conclusion, the work of Kojima et al 1995a confirms the theoretical predictions 
(Rice et al 1993) that a 3-leg system becomes magnetic at low temperature but not a 
2-leg system. 
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Dalmas de Reotier and Yaouanc: fig. 8b 



°- 0.5 




2 4 6 

Time (jxs) 



1.0 



Q- 0.5 



— i — i — i — | — i — i — i — | — i — i — i — | — i — i — 

Sr 2 Cu 4 6 (2-leg ladder) 

20 mK 



0.0 




2 4 6 

Time (us) 



Figure 8. Some spectra recorded on Sr2Cu406 which has a 2-leg spin ladder 
structure. The solid lines are fits (adapted from Kojima et al 1995a). Note that 
the horizontal scales are ~ 10 times larger than in figure |^ 



3.3. Diamagnetic domains in beryllium 

In general terms, we first note that the relation between induction, external field and 
magnetization for a given material is given by 

B = B oxt -,i (N-l)M, (16) 

where N is the demagnetization factor tensor and 1 the unit tensor. Therefore if 
the external field B cxt is applied perpendicular to an infinite plane, in practice to 
a plate like sample with extremely small thickness relative to radius, (16) yields (see 



section 2^) B = B ox t, i.e. the induction should follow the external field. 

It is well known that the magnetic response of a non-magnetic metal can oscillate 
as a function of i? cx t ■ This is the de Haas- van Alphen (dHvA) effect which is used to 
study Fermi surfaces (Shoenberg 1984). Condon 1966 noticed that if the oscillating 
amplitude is large within some part of each dHvA oscillation period, i.e. for a 
given -B ex t range, the conduction electron states in this range are thermodynamically 
unstable and cannot follow the electrodynamics relation B = B cxt . Therefore the 
electronic system should jump periodically over the forbidden intervals of B. A 
compromise could be the splitting of the magnetic energy in alternating diamagnetic 
(with induction smaller than B) and paramagnetic (with induction larger than B) 
domains. With such a domain structure the relation B = B cx t is fulfilled as an average 
over the sample. The so-called Condon domains are spectacular manisfestations of the 
collective behaviour of the electrons in quantized cyclotron orbits, i.e. Landau states. 

The first direct observation of the Condon domain formation was made on silver 
by NMR measurements (Condon and Walstedt 1968). One had to wait 28 years for 
a second report on the domain formation: using transverse field fjSR measurements, 
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Solt et al 1996a have reported the observation of domains in beryllium. Interestingly, 
the prediction of Condon was made for beryllium. 

As a first step, Solt et al 1996a have analyzed their data using (7). The field 
dependence of Xx is presented in figure [| Any deviation from Xx(B ex t) ~ constant 
must reflect the influence of the dHvA effect. The function Ax(-Boxt) rises periodically, 
reaching values about 10 times as large as its minimum. The period is consistent 
with the expected dHvA value. The broadening model is practical to visualize the 
dHvA effect but not strictly justified. Solt et al 1996a have then performed a Fourier 
analysis. The results for the central peak region of figure ^ are presented in figure |l(| 
As expected, within a given field range, B\ oc , which is a measure of the induction, 
does not follow -B cx t- The diamagnetic and paramagnetic signals are clearly seen, their 
respective populations changing smoothly as -B cx t increases. Interestingly, the domains 
could not be visualized by NMR on beryllium because of the excessive linewidth 
produced by the electric field gradients. Since the muon does not have a quadrupolar 
moment, it does not couple to these gradients. The Condon domains have been 
detected by a point-like probe such as the muon because the volume occupied by the 
domains is much larger that the volume of the domain walls: at -Bcxt = 1 T, the 
domain and wall thicknesses are ~ 30 /im and ~ 1 /mi, respectively (Solt et al 1996b). 
Note that the observation of the domains means that they are pinned (static) for the 
time scale of the muon experiment (~ 1 /is). 
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Figure 9. Exponential damping rate Ax of the /xSR signal as a function of the 
intensity of the external field measured on beryllium. The temperature is 0.8 K. 
The periodic sharp rises of Ax mark the onset of line splitting due to the domain 
formation. For Ax > 0.4 MHz the broadened line turns out to be a well resolved 
doublet. The solid line is a best fit to a truncated Fourier series (adapted from Solt 
et al 1996a). 



The Condon domains discovered in silver and beryllium correspond to the simplest 
possibility. A variety of domain-like periodic structures have been predicted (for a 
review, see Solt et al 1996b). We expect that with the improvement of the high field 
transverse field /iSR spectrometers, the investigation of magnetic domains in the bulk 
of non-magnetic metals will attract much interest in the future. 
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Figure 10. Analysis of the central peak of figure [J The (a) panel displays the 
intensity of the local fields at the muon site versus the intensity of the external field. 
The difference between the fields in the paramagnetic and diamagnetic domains is 
clearly resolved. For reference the line whose equation is Bi oc = B ex t is plotted 
(dotted line). The (b) panel presents the relative populations of the paramagnetic 
and diamagnetic domains versus the intensity of the external field. The lines are 
guide to the eyes, (adapted from Solt et al i996a). 



3.4- Localized and itinerant f electrons in heavy fermion materials 

In recent years the analysis ol /iSR data recorded in some intermetallic compounds 
containing cerium or uranium atoms has suggested that two different substates of 
electrons of / character are present. In this section the available experimental evidence 
found recently is presented. We discuss in some detail results obtained for the 
heavy fermion superconductor UPd 2 Al 3 and the heavy fermion metal CeRu2Si2- We 
then mention the superconductor CeRu2 and give a short discussion of the magnetic 
behaviour of these / electron magnets, pointing out the complementary nature of the 
information obtained by /iSR and inelastic neutron scattering. 

UPd2Al3 belongs to the small family of heavy fermion superconductors; see for 
example Heffner and Norman 1996. It has a hexagonal crystal structure (space 
group P6/mmm) and exhibits a coexistence of superconductivity (T c ~ 1.5 K) 
and antiferromagnetism (TV — 14 K). The uranium atom carries a relatively large 
magnetic moment of 0.85 /is as determined by neutron diffraction. The uranium 
moments are aligned ferromagnetically in the basal plane along the a axis and stacked 
antifcrromagnetically along the c axis. The magnetic structure is not affected by 
the superconducting transition, which is intriguing in view of the unquestionable 5/ 
character of both the heavy quasiparticles forming the Cooper pairs and the electrons 
responsible for the antiferromagnetic structure. The transverse held /iSR investigation 
of the superconducting phase by Feyerherm et al 1994 provides information on the 5/ 
electrons involved in the Cooper pair formation. 

Since the local field produced by the magnetic sublattices in the antiferromagnetic 
phase cancels at the muon localization site, a study of the muon frequency shift could 
be undertaken below T/v and in particular in the superconducting phase. The key 
point is the observation that the positive frequency shift along the c axis as well as 
the negative frequency shift in the basal plane are decreasing in absolute values when 
the temperature decreases below T c as seen in figure [ill In discussing the possible 
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origins for this decrease, Feyerherm et al 1994 noticed that the diamagnetic shift due 
to the flux expulsion or a possible change of the hyperfine coupling cannot explain 
the data. Therefore, the observed partial reduction of the frequency shifts can only 
reflect a decrease of the susceptibility of the 5/ electrons. Taking into account the 
strong anisotropy of the bulk susceptibility, these authors infer that the / electron 
susceptibility reduction due to superconductivity is isotropic. This, apparently 
surprising, result can be undertood if the 5/ electrons are viewed as two essentially 
independent electron subsets. The decrease of the susceptibility is associated with 
the electron system formed by the heavy quasiparticles condensing into Cooper pairs 
below T c . Interestingly, this reduction suggests singlet pairing of the Cooper pairs. 
The residual susceptibility is ascribed to the electron subsystem associated with the 
local antiferromagnetism which is not affected by superconductivity. We mention 
that the analysis of NMR Knight shift data leads to the same physical picture if one 
assumes that the hyperfine coupling constant is temperature independent, in contrast 
to the fj.SH analysis which does not need this assumption (Kyogaku et al 1993). 
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Figure 11. Temperature dependence of the muon frequency shift in UPd2Al3 
measured in an external field of 1 T applied either along or perpendicular to the 
c axis. Note the strong reduction of the absolute value of the shifts below the 
superconducting transition (T c ~ 1.5 K). The solid lines are guides to the eyes. 
The data are not corrected for demagnetization and Lorentz fields, (adapted from 
Feyerherm et al 1994). 

The picture of an isotropic itinerant 5/ electron subset is supported by the fiSK 
observation of an almost isotropic London penetration depth (Feyerherm et al 1994). 
It is noteworthy that specific heat measurements under pressure indicate that the 
itinerant subset accounts for 80 % of the linear coefficient of the specific heat (Caspary 
et al 1993). 

CeRu2Si2 has attracted considerable interest because it lies at the borderline of 
a magnetic instability, between long-range magnetic order and paramagnetic ground 
state (for a review, see Kambe et al 1996). 

In figure |l2] is presented the temperature dependence of the zero-field exponential 
damping rate measured on a CeRu2Si2 crystal (Amato et al 1993). The rate increases 
significantly below ~ 2 K, corresponding to an enhancement of the field spread at 
the muon site of the order of 20 /xT. The static nature of the local field is proved by 
the strong reduction of the damping rate measured in a longitudinal field as shown 
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Figure 12. Temperature dependence of the damping rates measured with a single 
crystal of CeRu2Si2- On the (a) and (b) panels are presented the rates recorded with 
the longitudinal and transverse geometry, respectively. These measurements show 
that CeRu2Si2 exhibits, in addition to a static magnetic field distribution below ~ 
2 K, a dynamical magnetic component persisting up to ~ 150 K. The solid lines are 
guide to the eyes (adapted from Amato et al 1993) 



in figure [lj. The exponential shape of Pz (t) in zero-field is surprising as one would 
expect to observe a parabolic function for such small damping rates. Amato et al 1994 
argue that the enhancement of the field distribution originates from a static electronic 
magnetic moment carried by the Ce atoms. The Ce magnetic moment is estimated 
to be ~ 10~ 3 [1b- Because of this extremely small value, this magnetism is probably 
itinerant in nature. Since the value is found using a localized magnet model, there 
is quite a large uncertainty on the value of the magnetic moment deduced from the 
data. The small longitudinal damping rate which seems to exist at high magnetic 
field suggests the presence of fast electronic spin fluctuations. The existence of such 
fluctuations is confirmed by the transverse field measurements shown in figure |l2|. 
These fluctuations, which are observed up to ~ 150 K, have been related to the short- 
range dynamical magnetic correlations detected below ~ 60 K by inelastic neutron 
scattering (Regnault et al 1988). They involve a relatively large Ce magnetic moment 
(0.6 (is, see Amato et al 1994) and an extremely short correlation time of ~ 10~ 13 s 
(Amato et al 1993). 

Recently, zero-field and longitudinal field measurements on a polycrystalline 
sample of the superconductor CeRu2 (T c = 6.1 K) have shown that this compound 
presents a magnetic phase transition at ~ 40 K characterized by an extremely small 
magnetic moment of ~ 10~ 4 (jlb (Huxley et al 1996). Because of this extremely 
small value, the magnetism is itinerant in nature. We note that a high energy 
neutron scattering study of CeNi2, a compound which should have the same magnetic 
characteristics as CeRu2, has revealed a strong inelastic paramagnetic response of 
the Ce ion with a large Kondo temperature (Murani and Eccleston 1996). The 
combination of the /xSR results on CeRu2 and the neutron data on CeNi2 suggests 
again the existence of a two component 4/ electronic system: in this case, a 4/ weakly 
polarized itinerant electron system and a paramagnetic 4/ localized electron system. 

The picture of two different substates of electrons of / character which emerges 
from the analysis of the UPd2Al 3 , CeRii2Si2 and CeRu2 data was already suggested 
for UCu5 by Schenck et al 1990. Some years ago inelastic neutron scattering 
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measurements on CeCug and CeRu 2 Si2 were analyzed in terms of two magnetic 
contributions (Aeppli et al 1986 and Regnault et al 1988). Amato et al 1993 have 
shown that the two contributions observed by neutron and /^SR can be nicely related 
in the case of CeRu2Si2- Clearly, it would be of interest to carry on this type 
of comparison keeping in mind that, contrary to the inelastic neutron scattering 
technique, the /xSR method probes the spectral weight of the modes at extremely 
small energy transfer, i.e., quasi-static modes. 

The two component picture which is proposed from the analysis of the data 
is in qualitative agreement with the duality model for heavy fermion (Kuramoto 
and Miyake 1990) which introduces two coupled electron components: itinerant and 
localized. The physical basis of such a picture is attributed to the fact that the one- 
particle density of states has a triple peak structure in strongly correlated fermion 
systems. This structure consists of two broad peaks corresponding to the upper and 
lower Hubbard bands and a narrow quasiparticle peak at the Fermi level. Interestingly, 
this type of structure is found by recent numerical works using the d — oo technique, 
d being the space dimension (see for example, Georges and Kotliar 1992). A recent 
theoretical discussion is given by Pepin and Lavagna 1997. 

A high resolution photoemission study of CeRu2 has found a substantial 4/ 
electron density at the Fermi level, in agreement with the interpretation of the fiSR 
result on this compound and the duality picture (Yang et al 1996). 



4. Spin dynamics in magnets 



The study of the spin dynamics in magnets has been a traditional subject of //SR. 
But it is only recently that quantitative information has been extracted, thanks to 
greatly improved experimental conditions and better data analysis. In section 4.1 and 
section |4.2j , spectra recorded respectively in the critical regime of a ferromagnet and for 
a spin-glass are analyzed in terms of spin-spin correlation-functions. The first study 
has taken advantage of the possibility to carry out the measurements in truely zero 
field and the second one, of the large spectrum of fluctuations which can be probed. 



The discussion is partly based on the material presented in Appendix B and Appendix 



4-1- Critical and low temperature spin dynamics in ferromagnets 

The first detailed analysis of data recorded in the critical regime has been given for 
Ni and Fe by Yaouanc et al 1993a and 1993b and for Gd by Dalmas de Reotier and 
Yaouanc 1994. These works were restricted to the critical paramagnetic regime. In 
this section we present an analysis of the spin dynamics of the Gd 3+ ion spins in the 
dipolar axial ferromagnet GdNi 5 covering the whole temperature range (Yaouanc et al 
1996a). This zero-field study reveals the effect of the dipolar interaction on the physical 
properties of the spin-spin correlation-tensor A"^(q), both in the paramagnetic and 
ferromagnetic state. This effect is expected to be strong since the measurements probe 
the fluctuation modes in the zone center where the long range nature of the dipolar 
interaction dominates the dynamics. This interaction introduces non-conserving terms 
in the Hamiltonian which prevents the slowing down of the longitudinal (to the 
wave vector) modes near the Curie temperature, Tq (Kotzler 1986). Therefore, if 
the measured relaxation rate reflects only these modes, it should saturate as the 
temperature approaches Tc- 
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GdNi 5 crystallizes in the hexagonal CaCus crystal structure and exhibits a 
ferromagnetic phase transition at Tq — 32 K characterized by a small magnetic 
dipolar anisotropy field as determined by magnetization measurements: B a (T = 5 
K) ~ 0.21 T. The muon localization site(s) is unknown. 

An overview of the zero-field relaxation rate \z is presented in figure |l3|a. This rate 
reflects the fluctuations of the Gd 3+ ion spins. No signal is found in the ferromagnetic 
state when is perpendicular to the crystal c axis since the spontaneous muon spin 
rotation is too fast to be resolved at the pulsed source (ISIS) where the data were 
recorded. Four temperature regions can be distinguished: far above or below the 
Curie temperature Tq and the critical paramagnetic and ferromagnetic regions. 

Far above T@, the wave vector dependent spin-spin correlation-tensor, A Q/3 (q), 
is expected to be isotropic and independent of the wave vector: A Q/3 (q) = K8 a P . 
The observed anisotropy of \z, which reflects the anisotropy of the coupling tensor 
between the muon spin and the Gd 3+ spins, can be explained if the muon occupies 
the interstitial site of coordinates (1/2,0,0). 



Dalmas de Reotier and Yaouanc: fig. 1 3a Dalmas de Reotier and Yaouanc: fig. 1 3b 
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Figure 13. (a): an overview of the zero-field relaxation rate of GdNis (Yaouanc et al 
1996a). The two sets of data refer to measurements performed on two crystals which 
differ by the orientation of S M relative to the c axis, (b): zero-field relaxation rate of 
GdNis at low temperature (Yaouanc et al 1996b). The full line is the prediction for 
the relaxation induced by the Raman spin wave process (equation (17)). 

Far below Tc, the spin-lattice relaxation rate increases with temperature as shown 
in figure [l3[b. Its thermal variation is interpreted as due to the relaxation of the 
muon spins by the magnons. Since the minimum magnon energy is much larger 
than the Zeeman muon energy htu^, a one magnon process cannot flip the muon 
spin because energy-conservation requirement would not be ensured. Therefore, the 
perturbation operator which induces the muon spin flip cannot contain the J x and 
J y operators (J denotes the Gd 3+ spin). In terms of the correlation-tensor A Q/3 (q), 
this means that one does not have to consider the correlations with {a(3} — {x,y}. 
On the other hand, a two-magnon, or Raman, process does not present problems 
as regards energy conservation because the only requirement is that the energies of 
the annihilated and created magnons must be equal (we neglect Tioj^). Therefore the 
energy principle tells us that only the parallel (to the easy axis z ; the Z and z axes 
are parallel) fluctuations contributes to the relaxation, i.e. the measurements only 
probe the correlation-function A 2Z (q) (Yaouanc and Dalmas de Reotier 1991). 
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While in NMR the Raman process is almost never observed because the hyperfinc 
interaction between the spin probe and the lattice spins is isotropic (therefore a 
perturbation operator such as a±J z does not exist; a is the Pauli operator of the 
spin probe), this process is expected to be active in /zSR since the interaction between 
the spins is mainly due to the dipolar interaction which is spatially anisotropic. 

The key parameter which determines if the muon spin is relaxed by the Raman 
process is the ratio of the minimum magnon energy over the thermal energy. The 
relaxation is effective only if this ratio is sufficiently small. This means that Xz = 
at very low temperature as observed experimentally in figure |l3|b. For a small ratio, 
i.e. at high temperature, we expect Xz oc T 2 since two magnons are involved. This 
quadratic law is a robust result in the sense that it does not depend on the details of 
the model. 

Using the simple magnon dispersion relation hu>q = D m q 2 + A, where D m is the 
magnon stiffness constant and A the anisotropy energy which is of dipolar origin, and 
the fact that A = gLHBB a <C ksT for GdNi5, the following approximate result holds 
(Dalmas de Reotier and Yaouanc 1995): 



\ CglT 2 
Az = 



D 3 



1 + y [(C-(q = 0)) 2 + (C^(q = 0)) 2 ] 



ln(fc B T/A). (17) 



C is the analytical part of the tensor describing the dipolar interaction between the 
muon spin and the Gd 3+ spins; see (A10). Equation (17) accounts for the effect of the 
spin waves near the zone center. Note that Xz is independent of the characteristics 
of the muon localization site(s) if C zx (q = 0) = C zy (q = 0) = as it is the case 
for the possible muon sites in GdNis. g^ is the Lande factor and Ub the Boltzmann 
constant. C is a universal constant (C = 129.39 (meV) 3 .A 6 .s _1 .K" 2 ). The fit to (17) 
presented in figure |l3| b yields D m = 3.2 (1) meV.A 2 . This estimate allows one to 
extract the value of an important parameter: namely the dipolar wave vector, qn, 
which determines the relative strength of the dipolar and exchange interactions. The 
analysis gives qn = 0.19 A -1 . The simple linear magnon theory cannot describe the 
data above ~ 20 K because the magnon-magnon interactions are neglected. 

We now analyze the paramagnetic critical behaviour of Xz- As shown in figure |l4|a, 
Xz is almost isotropic. This is expected for a dipolar Heisenberg paramagnet for which 
the correlation-tensor is 

A^(q) = A T (g)P^(q) + A £ (g)if*(q) (18) 

where P^ 7 (q) and P^ 7 (q) are the transverse and longitudinal (relative to q) projector 
operators, respectively. This decomposition reflects the symmetry of the Hamiltonian 
at small q value, i.e. of the long range dipolar interaction between the Gd 3+ ions 
(Frey and Schwabl 1988, 1989 and 1994). Since the dynamics near Tc is driven by 
the modes at small q, we need to consider the expansion of the muon-lattice coupling 



H is 



tensor G Q/3 (q) only near q = 0: G Q/3 (q) = -4tt P^(q) - C Q/3 (0) - SIM 
the hyperfine tensor. These functional forms of A and G, together with the expression 



for Xz given in Appendix C, leads to the simple result 



X z = W{a L I L (v)+a T I T ( V )}, (19) 

where W is a nonuniversal constant, cll,t depends only on the muon localization 
site(s) and I L ' T (<p) are universal fluctuation functions of the temperature through the 
angle ip. I L {(fi) and I T (<p) account for the longitudinal and transverse fluctuations, 
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respectively. We have <p = arctan(q£i£) with £ = £ > Qt~ v . £ is the correlation length, £ 
the correlation length extrapollated to T = 2Tq, t = \T—Tc\/Tc and v the correlation 
length critical exponent (y ~ 0.69). Therefore, in general, Xz is a weighted sum of 
I L {<-p) and I T {ip). Experimentally, a saturation is observed when approaching Tc (see 
figure [14ja). This is understood if cilI l {(p) ^S> ar^if) since I L ((fi) saturates as Tc is 
approached. 



Dalmas de Reotier and Yaouanc: fig. 14a Dalmas de Reotier and Yaouanc: fig. 14b 
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Figure 14. (a): zero-field relaxation rate measured in the critical paramagnetic 
state given as a function of the temperature relative to the Curie temperature and 
orientation of S M relative to the c axis (same symbol convention as in figure 
). The full and dashed lines are predictions of the mode coupling theory for the 
critical behaviour of Xz (T) in a dipolar Heisenberg ferromagnet (Frey and Schwabl 
1988, 1989 and 1994). (b): zero-field relaxation rate measured in the ferromagnetic 
state near the Curie temperature. The full line is the prediction for the critical 
paramagnetic fluctuations. The relative weight of the longitudinal and transverse 
fluctuations is taken as given by the analysis of the paramagnetic fluctuations. Both 
figures are from Yaouanc et al 1996a. 



Dalmas de Reotier and Yaouanc: fig. 15 



2.0 
1.5 


i i I i i I i i I i i i 






D 

^ 1.0 
cr 


: \ I 

s \ 


1 /a? 


0.5 






0.0 


: fenomagneticX 


/ paramagnetic 


i . , , , i . , , , i . , , , i . , , . i 





2 10 12 



Figure 15. Halperin-Hohenberg diagram for the /^SR measurements of GdNis at 
|7c - Tj/T c = 10 -2 (from Yaouanc et al 1996b). Their locations are indicated by 
the two bars. They have been clearly performed in the critical regime delimited by 
the shaded region. 
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An interesting result of this study is the observed similarity between the 
paramagnetic and ferromagnetic longitudinal critical fluctuations as discovered when 
comparing figure |l4^ , and figure |l4|b. Using the dynamical scaling theory of Halperin 
and Hohenberg (Halperin and Hohenberg 1967), it can be understood as follows. The 
basic quantity which distinguishes the different regions in the plot is the 

product q£. Since the measurements are mostly sensitive to longitudinal fluctuations, 
they probe the modes with q ~ qn (Dalmas de Reotier et al 1994). Therefore the 
relevant quantity is q^S,- Taking £ = 1 A one finds qo£ — 5 at t — 10~ 2 . Despite 
this rough estimate for the correlation length, the measurements are still clearly in the 
critical regime of the paramagnetic and ferromagnetic dynamics as shown in figure |l5|. 
Thus the continuity of the dynamical behavior when crossing Tq and therefore the 
observed similarity is understood. Nevertheless this argument calls for a detailed 
theoretical justification : in a dipolar magnet two scaling variables are needed instead 
of one for the isotropic model of Halperin and Hohenberg. 

We note that the relaxation rate of the points very close to Tc (namely the points 
which have been used for the determination of Tc and correspond to t < 0.0015) is 
significatively larger than the saturation value obtained for 0.002 < t < 0.02. This 
can be seen by comparing the values of the damping rate in figure |l4| a and figure |l4| b 
with figure [l^a. This increase of the damping rate near Tc could be due to the Ising 
crossover that has been observed in metallic Gd (Dalmas de Reotier and Yaouanc 
1994). 

We point out that systematic zero-field investigations of the dynamics in rare- 
earth metals would be of great interest: it could provide essential information for the 
determination of their universality class through the measurement of the dynamical 
critical exponent. 

Recently Heffner et al 1996 have measured Xz in the ferromagnet 
Lao.67Cao.33Mn03 which has attracted much interest because of its colossal mag- 
netoresistance. Equation (17) fails to predict the order of magnitude of Xz- This 
may not be surprising since this oxide is an inhomogeneous mixed-valence compound. 
The results of Heffner et al 1996 reveal the presence of a density of states for mag- 
netic excitations larger than expected, which has a strong influence on the correlations 
along the easy axis. Its origin could be related to the existence of a second d electron 
component in addition to the localized d component responsible for the conventional 
spin waves detected by inelastic neutron scattering (Perring et al 1996, Moussa et 
al 1996). This picture is identical to the one used to describe re sult s obtained on 
some metallic compounds containing Ce or U atoms (see section 3.4). In order to 
characterize these unexpected excitations further, it would be of interest to extend 
the measurements to crystals and investigate the possible field dependence of Xz over 
the whole temperature range. 



4-2. The correlation-function in spin-glasses 

The intensive investigation of glass and spin-glass forming systems in the eighties has 
shown that the key to a comprehensive understanding of the transition lies in the 
dynamics (Fisher and Hertz 1991). The spin-spin self correlation- function is the most 
important quantity in the spin glass systems since cross correlations are zero, i.e. one 
can neglect the wave vector dependence of the correlation-function. Therefore for an 
isotropic spin glass system one should only consider the correlation-function A(t) = 
A(t, q = 0) (see Appendix C). 
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Below the spin glass transition temperature T g , a fj.SH study has shown that 
A(t) decays as a power law after some microscopic time, of the order of 10~ 14 s 
(MacLaughlin et al 1983). However it is only recently that the form of A(t) above T g 
has been established (Keren et al 1996b). In this section we highlight this work. 

The goal of the work by Keren et al 1996b was to distinguish, using the longitudinal 
field method, between three possible forms of A(t) above T g : the power law (d- \t\~ a ), 
the stretched exponential (c?exp[— (C|£|)' 3 ]) and the cutoff power law (d ■ |£|~ Q /(CI*I)) 
which is often approximated by the Ogielski form (d ■ \t\~ a exp[— (CKI)' 3 ])- The 
exponents a and [3 are positive by definition. The difference between these forms 
is fundamental: the power law is time-invariant, the stretched exponential has a well- 
defined time scale given by l/£ and the cutoff power law is time-invariant only at 
times shorter than l/£. 

Since the longitudinal field method does not probe directly A(t) but only Pz(t), 
one has to determine the relation between these two quantities. Because the magnetic 
impurities are randomly distributed, the spin environment of each muon is different. 
We first consider an expression of Pz(t) for a given environment. As written in (10), 
for a crystal in the fast-fluctuation limit, Pz{t) is an exponential function with a 
relaxation rate \z = 2 A 2 /v. The field dependence of 1/v is given in terms of the 
correlation- function: 

1 r°° 

A(t) cos (7 M #oxti) dt. (20) 



i/(Boxt) 2A(0) J_ 



This formula, valid for an isotropic system, can be deduced from Appendix B. A(t) 
does not depend on B cxt if the electronic Zeeman energy is much smaller than the 
spin-spin coupling energy. The parameters A and v can vary from one muon site 
to another, hence their average values should be considered. Most treatments have 
assumed that v is site independent. In order to understand their data, Keren et al 
1996b are forced to use a weaker assumption. They write v(B ex t) — c l(B ex t) where 
the site dependence enters through the prefactor c. Thus the measured depolarization, 
which is an average, is given by 



P z (t) = / / p(A,c)exp 



2A 2 ' 
t 



c l(B cxt ) 



dc dA, (21) 



where p(A,c) is the probability of occurence of the values A and c. For the three 
possible forms of A(i), Pz(t) obeys asymptotically the scaling relation 

P z (t,B ext )=Pz(t/B2 xt ), (22) 

where 7 = 1 — a for the power law and the cutoff power law, and 7 = 1 + /? for 
the stretched exponential form. Asymptotically means for the stretched exponential, 
7/i-B ex t ^> C; an d for the cutoff form, <C 1. 

It is interesting to notice that (22) should be valid with 7 =1 for a conventional 
magnet if the spin dynamics is slow because A(t) is an exponential function for such 
a magnet. 

The measurements have been done on the canonical Heisenberg spin glass system 
AgMn(0.5 at. %). In figure [l6| a we present the field dependence of the spectra 
recorded at 3.2 K, i.e. just above T g . In figure |l6| b the validity of the scaling law 
(22) is demonstrated for 7 = 0.76 (5) over three orders of magnitude in t/B2 xt - The 
stretched exponential form is inconsistent with the measured 7 value. According to 
(22), an instantaneous relaxation should occur as B ext — > since then Pz(t) — » at 
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any t. This is not observed as confirmed in figure [i6|a. But as B cyit — > the exponential 
term of the cutoff power law provides a cutoff. Therefore, only the cutoff power law 
describes the field dependence of Pz{t)- This conclusion is achieved without assuming 
a specific functional form of Pz (t) . 



Dalmas de Reotier and Yaouanc: fig. 1 6a Dalmas de Reotier and Yaouanc: fig. 1 6b 
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Figure 16. (a) Dependence of the fj.SK spectra on the intensity of the longitudinal 
magnetic field recorded on the spin glass AgMn(0.5 at. %). The measurements have 
been performed at T = 3.2 K, i.e. above the spin glass transition temperature (T g = 
2.95 K). The solid lines are guides for the eyes, (b) The same spectra as presented 
on figure |l6p but plotted as a function of the scaling variable t/B^ t e and for various 
values of -B cx t. t is the time and -B e xt the external magnetic field. Both figures are 
adapted from Keren et al 1996b. 

A complementary approach is to test a functional form for Pz (t) . One possibility, 
which has been tested with success at a high impurity concentration (Campbell et al 
1994), is the stretched exponential function: 

Pzit) =exp[-(At)' 9 ] . (23) 

Keren et al 1996b show that this function provides a very good description of the 
AgMn(0.5 at. %) spectra recorded at various temperatures around the spin glass 
transition temperature by cooling the sample in a longitudinal field of 12.5 mT. In 
figure [n] are plotted the two parameters characterizing the stretched exponential. The 
maximum of A indicates the vicinity of the phase transition. As the temperature is 
lowered towards T g , the exponent f3 saturates to 1/3. This result holds for a wide 
range of concentration in AgMn. 

Keren et al 1996b have been able to distinguish between three possible forms 
of the spin-spin self correlation- function with only a scaling argument. Using the 
neutron spin echo technique, Mezei and Murani 1979 could not achieve this result. 
Interestingly, the success of the /xSR work lies in the use of the spin correlation- 
function concept for the data analysis. It would be of interest to check if the limit 
/3 — > 1/3 as T — > T g is a universal characteristics of spin glasses. 

5. Probing the magnetic properties of superconductors 

/iSR spectroscopy has been intensively used in recent years to probe the magnetic 
properties of superconductors. Below we present three experimental works which 
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Figure 17. Temperature dependence of the two physical parameters characterizing 
the stretched exponential used to fit spectra recorded on AgMn(0.5 a. %) at various 
temperatures around the spin glass transition temperature (T 9 = 2.95 K) by cooling 
the sample in a longitudinal field of 12.5 mT (adapted from Keren et al 1996b). The 
dashed lines are guides for the eyes. 



have greatly benefited from the availability of high quality crystals. In these works 
the physical properties of the mixed state are probed but with quite different results. 



In section |5.l| and section 5.2 the field distribution due to the vortices is visualized for 
YBa2Cu3 6.95 and Bi2+ x Sr2- x CaCu208+<5, respectively. Whereas the observation of 
a conventional 3D flux-line lattice for YBa2Cu306.95 leads to the analysis of the data in 
terms of the symmetry properties of the superconducting order-parameter, the strong 
temperature and field dependence of the Bi 2 + x Sr2_ :l: CaCu208+«5 field distribution 
offers the possibility of investigating its phase diagram in the temperature-field plane. 
In section |5~3| , the measurements in UPt3 show that, in addition to the analysis of the 
field distribution in terms of the symmetry of its superconducting order-parameter, 
one is able to unravel some basic magnetic properties of the compound and their 
interplay with superconductivity. 

5.1. The symmetry of the superconducting order-parameter in YBa2 CU3 06.95 

The symmetry of the superconducting order-parameter of the cuprate superconductors 
is the subject of ongoing reseach. A possible way of obtaining information on this 
symmetry is to investigate the excitation spectrum from the temperature dependence 
of the superconducting condensate density, n s . The flux line lattice in YBa2Cu3C>6.95 
is expected to be conventional, i.e. not to melt (at least for a usual magnetic field 



intensity; see section 5.2). Therefore the field distribution due to the flux line lattice 
is characterized by the coherence length and the London penetration depth, A, which 
is directly related to n s : namely 1/A 2 cx n s . An effective method to measure the field 



distribution is the muon spin rotation technique (see Appendix D). It probes the bulk 
of the material. 

Early /iSR measurements of the field distribution in sintered powders and low 
quality crystals concluded that 1/A 2 has a weak temperature dependence for T <C 
T c , suggesting the existence of an energy gap in the spectrum of excitations, in 
contradiction to the interpretation of other measurements. In particular, microwave 
data of Hardy et al 1993 on high quality crystals show that 1/A 2 b depends linearly on 
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1/2 

temperature. X a b is defined as (A a At) where A a and At, are the penetration depths 
for currents flowing along the a and b axis, respectively. This is consistent with the 
expected line of nodes for the order-parameter from singlet d x 2_ y 2 wave pairing. 

Although the microwave method has a high precision, it is not sensitive to the 
absolute value of A a fc(0) and probes only the skin depth at the surface. Sonier et 
al 1994 and 1997 and Riseman et al 1995 have performed /iSR measurements in a 
mosaic of single crystals. In figure [l^ we present the Fourier transformation of two 
spectra which have been recorded using a field cooling procedure with B cxt — 0.500 T. 
But while the distribution on the left has been recorded with that field value, the 
distribution on the right was taken in a field of 0.489 T after initially cooling the 
sample in a field of 0.500 T. Whereas the signal from muons in the sample remains 
practically unchanged due to the strong pinning of the vortex lattice, the sharp peak 
attributed to the background (muons stopped in the sample holder and cryostat walls) 
shifts down by ~ 11.2 mT. 



Dalmas de Reotier and Yaouanc: fig. 18a 
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Figure 18. Typical field distributions recorded in YBa2Cu306.95 after cooling the 
sample in a field of 0.500 T (adapted from Sonier et al 1994). While the distribution 
of figure |l8p has been recorded at 5.4 K with that field value, the distribution of figure 
p^ b was measured at 6 K after decreasing the field to 0.489 T. Because of the high 
superconducting temperature of this compound, the difference in temperature for 
the two measurements is expected not to be significant. The sharp lines below 0.50 
and 0.49 T in figure ^j| a and p^ b respectively represent the background contribution, 
i.e. muons which are not stopped in the sample. The other part of the distribution 
arises from the superconductor and is not affected by the field decrease. These 
measurements prove that the vortices are pinned. 



From the measured field distribution, Sonier et al 1994 and 1997 have extracted 
X a b using the theory explained in Appendix D, with the restriction that they do not 
use a proper cutoff function to account for the finite size of the vortex cores (Yaouanc 
et al 1997a). In figure Owe present 1/A^(T) for B cxt = 0.2 T, 1.0 T and 1.5 T. The 
linear dependence of 1/A^ 6 vs T, confirms the zero-field microwave measurements. We 
note that the weak field-dependence of l/A^ b (T = 0) can probably be explained using 
a proper cutoff function (Yaouanc et al 1997a). In their analysis, Sonier et al 1997 
take into account the in-plane anisotropy of the penetration depth discovered by fiSR 
(Tallon et al 1995, Bernhard et al 1995b) and infrared measurements (Basov et al 
1995). 

Measurements carried out by Riseman et al 1995 for 1.9 T < B cxt < 6.5 T show 
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Figure 19. Temperature dependence of 1/A^ b measured on YBa2Cu3 06.95 for 
three values of the external field applied along the c axis (adapted from Sonier 
et al 1997). The measured linear temperature dependence is consistent with a d- 
wave superconducting order-parameter. The observed field dependence is probably 
explained if a proper proper cutoff function is used (Yaouanc et al 1997a). 



that the Ginzburg-Landau paramater k is constant between 30 K and 75 K with a 
value k w 70. In addition, these authors finds an upper critical field B c2 = 90 (10) T. 

As pointed out by Sonier et al 1994, the sample quality seems to have a significant 
influence on the experimental results. This is understandable since A is in fact 
expressed in terms of an integrated excitation spectrum over parts of the Fermi surface 
(see Gross et al 1986 and Gross- Alltag et al 1991). Therefore, additional methods of 
investigation of the symmetry of the order-parameter are highly desirable. A recent 
example is given by Bernhard et al 1996 who have determined by ^iSR the depression 
of Tig cLS <X function of the level of Zn doping. They argue that the initial decrease 
of n s is inconsistent with s-wave pairing and magnetic scattering but points rather 
towards cZ-wave pairing, in agreement with the microwave data of Hardy et al 1993, 
the /xSR data of Sonier et al 1994 and 1997 and the micro-SQUID results of Tsuei 
et al 1994. However, this conclusion is disputed by Nachumi et al 1996. Therefore, 
untill the dispute is settled, it is not possible to give a definite conclusion concerning 
the interpretation of the effect of Zn doping on n s . 

Since the discovery of the high T c superconductors, trends in the relations between 
some of their parameters have been looked for. Some years ago a remarkable 
empirical relationship between T c and the low-temperature variance of the vortex field 
distribution was found for some high T c oxides (Uemura et al 1991). This experimental 
result was taken as evidence for a high-energy-scale pairing mechanism, consistent with 
the picture of real-space paired bosons. It was suggested that this scaling was valid 
for all the high T c materials. This suggestion has triggered an important experimental 
activity. The recent reports (Uemura et al 1993, Niedermayer et al 1993, Weber et 
al 1993, Bernhard et al 1995b, Zimmermann et al 1995, Tallon et al 1995) show that 
the simple universal scaling-law is only partially valid. For example, T c is not solely 
a function of the hole concentration: the oxides with and without Cu chains do not 
belong to the same class. 

We note that practically all these works have been performed in powder samples 
and the ^SR spectra have been analyzed supposing a Gaussian field distribution. 
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Although this methodology may not be completely safe (Sonier et al 1994 and 
Harshman and Fiory 1994), it has nevertheless been able to provide the first proof of 
the in-plane anisotropy of the London penetration depth in YBa2Cu30 2 _a (Tallon et 
al 1995). 

From the /iSR technical point of view, the works of Sonier et al 1994 and 1997 
and Riseman et al 1995 are remarkable since they have been performed both at high 
fields and with small samples (one of the Sonier's sample covers an area of 5 x 5 
mm 2 perpendicular to the c axis). In a near future one may foresee routine /iSR 
measurements on even smaller samples and therefore on samples of even higher quality. 



5.2. The vortex state in highly anisotropic high T c superconductors 

The state of the vortices in some of the highly anisotropic superconducting oxides is 
still a subject of discussion (see Bishop 1996a and 1996b, Nelson 1997 and Crabtree 
and Nelson 1997 for reviews). In these compounds such as Bi 2+a; Sr2_ x CaCu208+ < 5, 
the vortices are best described as layered systems of 2D pancake vortices interacting 
via a combination of tunneling Josephson currents and electromagnetic interactions. 
Increasing the field or the temperature, one expects to observe changes in the typical 
3D flux line lattice field distribution due to either disordering of the vortex lattice, a 
reduction of its dimensionality or its melting. 

In this section we focus on results obtained by Lee et al 1993, 1995 and 1997, 
Aegerter et al 1996 and Bernhard et al 1995a. 

At low temperature and for a field less than the crossover field B cr , a lattice of 



extended flux lines is observed. As seen in figure 18 and figure |2Cj, its /iSR signature is 
a typical strongly asymmetric field distribution with a pronounced tail towards high 
fields, showing that some of the muon spins precess in the local field caused by flux 
cores that are extended in the c direction. This interpretation is supported by the 
occurence of Bragg peaks in the neutron scattering experiments (Cubitt et al 1993). 
Drastic changes occurs when B cxt exceeds B cr . The field distribution becomes more 
symmetric as shown in figure Simultaneously the neutron Bragg peaks disappear 
indicating that the long-range coherence of the flux lattice is destroyed. The change 
of the shape of the field distribution is quantified by the skewness parameter a (see 



Appendix D| ). Its field dependence is presented in figure |2^. a is definitively smaller 



above B cr than below, reflecting the truncation of the high-field tail. This reduction 
is either due to the motion of the vortices (Harshman et al 1991, Inui and Harshman 
1993), a reduced dimensionality of the vortex structure (Brandt 1991 and Harshman 
et al 1993) or a transition to a glass phase (Ryu et al 1996, Ryu and Stroud 1996 
and Gingras and Huse 1996). The last two statements have common features and it 
may not be easy to distinguish between them (Gingras and Huse 1996). Lee et al 
1993 suggest that their data points to a reduction of the dimensionality of the vortex 
system: in a high field, the system consists of an array of pancake vortices that are 
uncorrclatcd in the c direction but ordered two-dimcnsionally within each stack of 
CuC>2 planes. This interpretation is supported by the numerical results of Schneider 
et al 1995. Such a dimensional crossover has been predicted to occur at B cr (Vinokur 
et al 1990 and Glazman and Koshelev 1991). Additional measurements (Bernhard et 
al 1995a) on underdoped and strongly overdoped crystals support this interpretation. 
Recently Aegerter et al 1996 have discovered that B CI = ^o/\ 2 ab where $o = 2.07 
x 10~ 15 Tm 2 is the quantum of flux. This latter result is understood in terms of a 
system of vortices controlled predominantly by electromagnetic interactions. This is 
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in contrast to other materials, such as Hgl^CaaC^Oio+a, where Josephson coupling 
plays a more significant role. 



Dalmas de Reotier and Yaouanc: fig. 20a 



Dalmas de Reotier and Yaouanc: fig. 20b 



1 1 i 1 1 1 1 i 1 1 1 1 i 1 1 1 
Bi 2 .isSr 1 . 8S CaCu 2 8 + 8 



(a) 




1.2 



CD 

I 1.0 



0.8 



1 i 1 1 1 1 I , 1 i J 1 1 i 1 1 1 i I , 1 i i 1 1 1 
Bi 215 Sr 185 CaCu 2 8 + 5 W 


• 

• 1 

.|B„ 

1 




1 

L 

1 \ 


• • 


r 




i • 

i 

i , , , , i , , , , 


• 

, , , i , 



50 100 150 200 
External field B ext (mT) 



250 



Figure 20. The vortex state in the highly anisotropic 

material Bi2.15Sr1.85CaCu2Og.1_15 as the intensity of the external magneti c fi eld is 
increased. The temperature is kept fixed far below T c . As shown in figure UOja, the 
shape of the field distribution becomes more symmetric at high field. The field de- 
pendence of the skewness parameter, which allows one to quantify the change in the 
shape of the distribution, is presented in figure pc| b (Aegerter 1997 and adapted from 
Aegerter et al 1996). Note the sharpness of the crossover at B CT . 

There is much interest in investigating the vortex state as a function of the 
temperature since it is expected that, at sufficiently high temperature, the vortex 
lattice should melt (Bishop 1996). In figure pl| we compare field distributions recorded 
at low and high temperature for a given field intensity. The high temperature 
distribution is very narrow and the shapes at low and high temperature are distinctly 
different. The temperature dependence of the skewness parameter indicates a sharp 
change at T m . a is even negative above T m , whereas in figure |2^ it is always positive. 
Comparing these results with the numerical calculations of the field distribution by 
Schneider et al 1995, T m is found to be the fusion temperature of the vortex lattice. 
A thorough analysis is presented by Lee et al 1997 which supports this interpretation. 
In addition, this analysis points out the determinant role of the electromagnetic 
coupling between the superconducting layers, In fact, it is this coupling rather than 
the Josephson coupling which determines the phase diagram below ~ 70 K, a rather 
high temperature relative to the superconducting temperature (T c ~ 80 K). 

The detailed characterization of the vortex state in the highly anisotropic oxides 
is a complex problem. The microscopic methods like /iSR and small angle neutron 
scattering, in combination with magnetization measurements, can lead to an improved 
understanding of the physics involved. For example, in Fi^S^CaC^Os+a, the "second 
peak" observed in the hysteresis loops and the sharp change in the field distribution 
detected by //SR. both occur at the same field, the crossover field B cr . allowing the 
origin of the "second peak" to be understood (Bernhard et al 1995a). 

5.3. Anisotropy of the magnetic response in UPt^ 

The hexagonal heavy fermion superconductor UPt3 has the unique physical property 
of showing two Meissner phases, at T c \ ~ 0.48 K and T c2 ~ 0.53 K, and three 
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Dalmas de Reotier and Yaouanc: fig. 21a 
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Dalmas de Reotier and Yaouanc: fig. 21b 
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Figure 21. The vortex state in the highly anisotropic mate- 

rial Bi2.i5Sri.85CaCu20g_|_^ as the temperature is increased. The intensity of the 
external magnetic field is kept fixed. As shown in figure Elk, the shape of the field 
distribution changes and becomes narrow at high temperature. The temperature 
dependence of the skewness parameter is presented in figure [2l| b (Aegerter 1997 and 
adapted from Lee et al 1997). Note the sharp change of a at T m and the negative 
value of a for T > T m . 



flux phases. In addition, neutron and magnetic X-ray diffractions have indicated 
an antiferromagnetic phase transition at T/v ~ 6 K characterized by a tiny magnetic 
moment of 0.02 (1) //s/U-atom at low temperature, lying in the basal plane along the 
b axis. This transition has never been detected by macroscopic measurements. The 
origin of the magnetic and superconducting phases of UPt 3 is one of the most debated 
subjects in condensed matter physics. The superconducting multiphases identify UPt3 
has a candidate for unconventional superconductivity. The term "unconventional" 
refers to the fact that the order-parameter has a lower rotational symmetry in the 
superconducting phases than in the normal state. More information on UPt3 can be 
found in the reviews of Sauls 1994 and of Heffner and Norman 1996. A review on 
unconventional superconductivity has been recently published (Muzikar 1997). 

We first present some recent zero-field measurements (Dalmas de Reotier et al 
1995). They have been performed with two purposes: to detect the magnetic phase 
transition at T/v and a signature of an eventual internal magnetic field induced by 
the Cooper pairs in the low temperature superconducting B phase (below T c i). The 
results are presented in figure |2^. For the large temperature range probed and the two 
crystals investigated (S^ either parallel or perpendicular to the c axis), Pz(t) is simply 
described by a Kubo-Toyabe function. The values of the damping rates are consistently 
explained as arising from the 195 Pt nuclear magnetic moments. Therefore neither the 
phase transition at T/v, nor an additional magnetic field in the superconducting B 
phase is observed. 

Since neutron diffraction measurements performed on some of the slices of the 
/iSR samples show that the antiferromagnetic phase transition still exists, one must 
conclude that either the electronic dipolar fields at the muon site cancel out exactly or 
the electronic uranium magnetic moments fluctuate too fast to be detected, i.e. their 
characteristic time is shorter than ~ 10~ 7 s. From neutron diffraction (Aeppli et al 
1988) we know that this time is longer than ~ 10~ n s. 

The zero-field measurements (Dalmas de Reotier et al 1995) indicate that a 
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Dalmas de Reotier and Yaouanc: fig. 22a 



Dalmas de Reotier and Yaouanc: fig. 22b 
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Figure 22. (a): Typical zero field spectra measured for UPt3. The full lines are 
fits. This figure shows that there is no additional local field at the magnetic phase 
transition and in the low temperature superconducting phase, (b): Temperature 
dependence of the Kubo-Toyabc damping rate A. The dashed straight line indicates 
its average. These results show that A is independent of the temperature and 
orientation of the crystal axes relative to S^. Both figures are from Dalmas de 
Reotier et al 1995. 



possible change in £>i oc induced by magnetism or superconductivity, if any, has 
to be smaller than approximately 3 /iT. Recently high precision magnetization 
measurements have shown that the possible change in bulk magnetization (rather 
than in the local field as observed by /xSR) in the superconducting B phase is smaller 
than 0.2 fiT (Kambara et al 1996). A Cooper pair produces an orbital magnetic field 
at the muon site. If it has a spin, a spin density has to be added to the orbital density. 
Since the orbital moment is expected to be much smaller than the spin moment, the 
measurements put a limit on the possible value of the spin density at the muon site. 
Taking into account the available theoretical estimates for the spin density, Dalmas 
de Reotier et al 1995 conclude that their results do not support models predicting a 
triplet spin state for the Cooper pair. 

Lussier et al 1996 have investigated the magnetic field response of UPt3 by single 
crystal neutron diffraction. Their results show that a field in the basal plane of up to 
3.2 T has no effect on the magnetic Bragg peaks. Since the intensity of these peaks is 
extremely small, the precision of the measurements is limited: it is only known that the 
angle of rotation of the magnetic moment in the hexagonal plane is smaller than 26°. 
Taking into account that the /iSR technique is well adapted to study small moment 
systems, Yaouanc et al 1997c have performed transverse high field ^iSR mesurements. 
They detect the magnetic phase transition if a large field is applied in the basal plane. 
As expected, the phase transition is not observed for a field applied along the c axis. 
Therefore the magnetic properties of UPt3 are found to be field dependent if probed 
at high field by a sufficiently sensitive method. 

The magnetic properties of UPt3 are remarkably sensitive to alloying. In the case 
of the pseudobinaries U(Pti_ x Pd x ) 3 , large uranium magnetic moments have been 
detected for 0.02 < x < 0.07. Neutron diffraction experiments on U(Pto.95Pdo.o5) 3 
show that the size of ordered magnetic moment equals 0.6 (2) /zs/U-atom, much larger 
than in the pure sample (0.02 (1) /zs/U-atom). In the zero-field /iSR measurements 
of deVisser et al 1997 on this compound, two frequencies (or one frequency and 
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a strongly damped Kubo-Toyabe signal; Amato 1997) are observed, indicating two 
magnetically inequivalent muon stopping sites. On the other hand, the zero-field 
signal in U(Pto.ggsPdo. 002)3 and in the pure sample are the same, i.e. no electronic 
magnetic signal is observed. 

The preliminary data of deVisser et al 1997 demonstrate a salient difference in 
the /iSR response between compounds with small and large magnetic moments in the 
U(Pti_ a; Pd a; ) 3 series. It is still not possible to extract reliable physical information 
from the limited amount of available experimental /iSR data. An important issue that 
could be resolved by combining these measurements with neutron measurements on 
the same samples concerns origin of the large difference between the intensity of the 
low temperature magnetic moment for UPt 3 and U(Pto.95Pdo.o5) 3 - 

The temperature dependence and the anisotropy of the magnetic field penetration 
lengths in the superconducting phases of UPt3 have been studied by muon spin 
rotation measurements (Yaouanc et al 1997b). The analysis of the temperature 
dependence of the penetration length parallel and perpendicular to the c axis has 
shown that the superconducting order-parameter in the B phase cannot just have 
a line of nodes in the equatorial plane of the Fermi surface. The analysis supports 
an hybrid order-parameter with point nodes at the poles and a line of nodes at the 
equatorial plane. A remarkable result of the measurements is the observation of a 
strong axial anisotropy of the penetration length. As explained at Appendix D|, t his 
anisotropy is directly related to the opening angle of the vortex lattice. In figure [23| we 
present the temperature dependence of this angle for B cxt applied along the a axis. 
If the penetration lengths were isotropic, it would be equal to 60°. It has already 
been measured by small angle neutron scattering (Kleiman et al 1992). But, because 
the neutron cross-section decreases dramatically as the temperature is increased, it 
could only be measured at 0.05 K. We note that the angle is pratically temperature 
independent. When B cxt is applied along the c axis of the crystal lattice, the angle is 
found to be temperature independent with a value of ~ 60°, i.e. the vortex lattice is 
a traditional hexagonal lattice. 
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1 1 1 1 1 1 1 

: + + + + + 
: UPt 3 


. - b - 


f B ext =18mT 
1 . .1 . .1 


1 , , , , 1 



0.0 0.1 0.2 0.3 0.4 
Temperature (K) 



Figure 23. Temperature dependence of the angle characterizing the vortex lattice 
of UPt3 for the external field applied along the a axis. The angle of the real space 
vortex lattice is defined in the insert, (adapted from Yaouanc et al 1997b). 
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6. Summary and perspectives 



In this review we have attempted to detail the possibilities offered by the /iSR 
technique for the study of the magnetic properties of materials. A quick look at 
all the examples described in this work shows that the prominently used properties of 
the technique are its high sensitivity to local fields, and its capability to characterize 
the time scale of their fluctuations. Interestingly, information can be extracted from 
compounds which do not display long range magnetic correlations. In fact, many /iSR 
investigations concern disordered or frustrated magnetic compounds. Two examples 
are given in this review. The typical spot size of a muon beam is some millimetres at 
a continuous source. Although this means that extremely small samples can not be 
studied, this macroscopic size allows one to probe magnetic structure with nanometre 
or micrometre scales. We have presented four examples of such studies. 

The fj.SH technique has many attractive feactures which explain its present 
fast development. Although we have not explained in any detail how a fiSK 
experiment is performed (this information is found in the general references provided 
at the beginning of section ^), relative to other non local microscopic techniques, a 
/iSR experiment is relatively straightforward. The /iSR spectroscopy is of general 
use: the sample does not need to contain special nuclei. A relatively involved 
sample environment can be used as illustrated, for example, by the high number 
of measurements performed in the 30 mK -IK range using dilution refrigerators. 
Measurements have been performed up to _B cx t = 6.5 T (Riseman et al 1995). The 
recording of a /iSR spectrum requires in most cases less than an hour. Therefore, in 
a reasonable amount of time, it is possible to map out the temperature and magnetic 
field ^tSR response of a compound. 

As does any experimental technique, the /iSR method has its own limitations. 
One of the major drawbacks is that the sample has to be sufficiently big. If a disk 
is used and if the muon beam is perpendicular to its plane, in the most favourable 
cases, the diameter of the disk should not be smaller that ~ 5 mm. A larger diameter 
is unfortunately required at a pulsed source such as ISIS, where the diameter has to 
be at least 20 mm. For the more widely used beamlines, the so-called surface muon 
beams (i.e. relatively slow muon beams), 150 mg of material per cm 2 of beam cross 
section is required to stop the muons. This means that the sample thickness needs not 
to be much larger than ~ 1 mm for organic materials (density > 1) and even much 
less (one to three hundred micrometres) for denser material such as oxides or metals. 

A priori, the muon diffusion and localization site(s) properties are unknown. This 
can turn out to be a problem to extract detailed information from the measurements. 
In this respect, it has recently been shown that frequency shift measurements may be 
of great help (see section 2.3). 

In order to respond to the demand of the users, major technical developments are 
scheduled at the ^SR facilities. 

The investigation of the vortex state in superconductors is at present one of the 
main subject of /iSR research. Optimized spectrometers in transverse geometry with 
intense magnetic fields at relatively high temperature are needed for this type of work. 
They should be available in a few years time both at TRIUMF and PSI. 

Surface science is a major topic of modern physics which, so far, is untouched by 
/iSR. The reason is that the momentum of the muons in presently available beam is 
too large and therefore the muons are stopped at best in a few hundred micrometres. 
Ultra slow muon beams are being developed at facilities. The challenge is to keep a 
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high degree of polarization in the slowing down process together with a reasonable 
efficiency. The possibility to moderately reaccelerate the ultra slow muons provides 
the opportunity of fine tuning the stopping range of the muons. For instance, such a 
beam will allow the study of magnetic multilayers which are of interest both for basic 
research and technological applications. 

The possibility of building an intense pulsed neutron facility in Europe, the 
European Spallation Source, is also being considered. As at the ISIS Facility, where 
pulsed neutrons and muons are available, many muon beams are scheduled to be 
introduced at this future facility. The experience gained at ISIS has shown that 
having beams of both particles at the same institution is a real advantage: it provides 
opportunities for the neutron and muon communities to interact. 

7. List of works published from July 1993 

The tables provide a list of published works. When a given work has been described 
in short articles and later on in an extended paper, we only refer to the latter report. 
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Table 1. A selection of work carried out on conventional magnets 



Chemical formula 



Comments 



References 



Magnetic correlations in magnets 



(Mm-^FesOPta 
MnSi 

Agi-^Mn^ 
Gd 

Ni and Fe 

Gdo.696Yo.304 
CeB 6 

PrCo 2 Si 2 
Y 6 (Mni_ a; Fe 2 ,)23 
CeSb 
EuO, EuS 

RbMnF 3 

Yo.97Sc .o3Mn 2 

Al-Mn-Si quasicrystal 

UNiGa 

YMn 2 

Ni 

Cri-^Fe^ 
/3-(NH 4 ) 2 FeF 5 
MnF 2 
CeRh 3 B 2 

Lao.67Cao.33Mn03 

YMn 2 Da; 

RcGae 

GdNi 5 

Ho 



dynamical crossover from ferromagnetism to 
asperomagnetism 

strong correlation between the 55 Mn electric 
field gradient and magnetic susceptibility 

spin correl atio n- function in a spin-glass 



1.2 



see section 

critical paramagnetic fluctuations and 
correlation-functions; effect of dipolar interaction 

evaluation of the Brillouin zone probed by /tSR 
study of critical paramagnetic fluctuations 

magnetic phase diagram investigation 

field response in disagreement with the 
hypothesis of an antiferroquadrupolar phase 

spin dynamics in an axial magnet 

static and dynamical magnetism in competition 

slow spin dynamics below Tjv 

theoretical prediction of the thermal behaviour 
of the relaxation rate; see Yaouanc ef al 1993b 

theoretical prediction of the thermal behaviour 
of the relaxation rate 

a spin liquid ground state is suggested 

study of the spin-glass state 

magnetic phase diagram investigation 

comparison between neutron and fiSK results 

temperature and pressure dependence of the 
Fermi contact field 

spin correlation-function in a spin-glass 

combined Mossbauer and fiSH analysis 

first proof of Ising dynamical critical behaviour 

investigation of the magnetic state 

ferromagnet with colossal magneto resip tance; 
unusual spin dynamics; see section 

effect of deuterium on spin dynamics 
and magnetic phase diagram 

combined fiSR and neutron study of these 
quasi-two-dimensional magnets 

effect of the dipolar interaction on the crit ical 
ferromagnetic fluctuations; see section 



1.1 



study of the incommensurate helicoidal structure 



Barsov et al 1993 and 1994 

Kadono et al 1993 

Campbell et al 1994 
Keren et al 1996b 

Dalmas de Reotier and 
Yaouanc 1994 

Dalmas de Reotier et al 1994 

Eccleston et al 1995 
Feyerherm et al 1995c 

Gubbens et al 1995 
Kilcoyne and Telling 1995 
Klauss et al 1995 
Lovesey and Engdahl 1995 

Lovesey et al 1995 

Mekata et al 1995 
Noakes et al 1995 
Prokes et al 1995 
Rainford et al 1995a 
Stammler et al 1995 

Telling and Cywinski 1995 
Attenborough et al 1996 
Brown et al 1996 
Cooke et al 1996 
Heffner et al 1996 

Latrochc et al 1996 

Lidstrom et al 1996a 

Yaouanc et al 1996a 

Krivosheev et al 1997 
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Table_2. A selection of work carried out on conventional magnets (continuation of 
Table |) 



Chemical formula 



Comments 



References 



Magnetic phase diagrams in organic magnets 



(TMTSF) 2 -X 

X = PF 6 , NO3, CIO4 

p-NPNN 

p-PYNN 

3-QNNN 

a-(BEDT-TTF) 2 KHg(SCN) 4 
1-NAPNN, 2-NAPNN 

MEM(TCNQ) 2 
p-CNPNN, 4-QNNN 



collective low-energy spin-density-wave 
in addition to single-particle excitations 

3D Heisenberg ferromagnet with Tq ~ 
effect of a longitudinal field; see section 



K; 



magnetic p 
see section 3.1 



lase transition at ~ 0.09 K 



0.21 K 



magnetic p hase transition at 
see section |3,l| 

detection of two spin density wave transitions 

1-NAPNN has a magnetic transition below 0.1 K 
and 2-NAPNN does not order magnetically 

spin dynamics in this spin-Peierls system 

p-CNPNN has a magnetic transition at ~ 0.17 K 
and 4-QNNN does not order magnetically 



Le et al 1993a 

Le et al 1993b 
Blundell et al 1995 

Blundell et al 1994 

Pattenden et al 1995 

Pratt et al 1995 
Blundell et al 1996 

Blundell et al 1997a 
Blundell et al 1997b 



Borocarbide materials 



YNi 2 B 2 C 


characterization of the vortex state of this 
conventional superconductor 


Cywinski et al 1994 


YNi 4 BC . 2 


proof that it is not a bulk superconductor 


Siillow et al 1994 


TmNi 2 B 2 C 


coexistence of magnetism and superconductivity 


Cooke et al 1995 
Le et al 1995 


ErNi 2 B 2 C 


observation of one frequency below Tjy 


Le et al 1995 


SmNi 2 B 2 C 
YNi 2 B 2 C 


exhibits a magnetic phase transition 
no magnetic correlations 


Prassidcs et al 1995 


HoNi 2 B 2 C 


observation of two magnetic transitions : one 
commensurate and one incommensurate 


Le et al 1996a 
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Table 3. A selection of work carried out on strongly correlated electronic systems 



Chemical formula 


Comments 




References 


Non-superconducting compounds 



CeRu2Si2 
Yi- x V x Pd y 

CeCus.gAuo.i 

CeCus.sAuo.s 
CeNiSn 



CeTrSn, Tr = Pt, Pd 

CePt 2 Sn 2 

UNi 4 B 
CeRhSb 
CeCus-ajAl* 
UCus-aPd* 

YbAuCu 4 

Cei-^Re^Nii-yTrySn 

CeTr 2 Sn 2 

Tr = Cu, Pt, Pd 

YbPdSb 

CePt 2 Si 2 

CeRuSi 2 



a two component 4/ system; see section 3.4 

competition between Kondo, RKKY and crystal field 
interactions: effect on the magnetic phase diagram 

no magnetic ordering in this non-Fermi-liquid 
system; Kondo disorder is maybe negligible 

phase transition at Tjv = 0.95 K 
magnetic correlations with surprising temperature 
dependence (the analysis yields a too small Kondo 
temperature; see Dalmas de Reotier et al 1996) 

CePdSn exhibits a simple second order transition; 
the transition in CePtSn is unusual 

investigation of magnetic correlations; see 
comments of Dalmas de Reotier et al 1996 

study of the magnetic phase diagram 

spin correlations below ss 0.6 K 

magnetic phase diagram investigation 

non-Fermi-liquid alloy; strong Kondo disorder 

effect of crystal field on Kondo-type fluctuations; 
possible muon induced crystal field effect 

effect of doping on the electonic ground state 

Investigation of the magnetic correlations in 
heavy fermion antiferromagnets 

a spin liquid system 

Kondo lattice compound with no magnetic ordering 

ferromagnetic phase transition at Tq = 11.6 K; 
very small ordered Ce magnetic moment 



Amato et al 1993 and 1994 
Wu et al 1994 

Amato et al 1995 
Bernal et al 1996 
MacLaughlin et al 1996 

Chattopadhyay et al 1995 

Kalvius et al 1995a 



Kalvius et al 1995b and 1995c 

Luke et al 1995 

Nieuwenhuys et al 1995 

Rainford et al 1995b 

Wiesinger et al 1995 

Bernal et al 1996 
Maclaughlin et al 1996 

Bonville et al 1996 

Flaschin et al 1996 
Lidstrom et al 1996b 

Bonville et al 1997 

Dalmas de Reotier et al 1997 

Duginov et al 1997 



Superconductors 



UPd 2 Al 3 

Cei +;E Cu2+j,Si 2 

UPt 3 

Ui_ I ,Re I Ru 2 Si 2 
Re = La, Y 

CeRu 2 

U(Pti_ I Pd a; )3 



isotropic; 



3,1 



London penetration depth approximate! 
a two component 5/ system; see section 

competition between superconductivity and magnetism 



absence of electronic magnetic signal in zero-fiek, 
study of London penetration depths; see section 

magnetic phase diagram investigation 



5.3 



detection of a magnetic phase transition with 
very small magnetic moments; see section 3.4 

study of effect of doping on magnetism; 
see section 5.3 



Feyerherm et al 1994 

Luke et al 1994 
Feyerherm et al 1995a 

Dalmas de Reotier et al 1995 
Yaouanc et al 1997b 

Cywinski et al 1995 
Park et al 1996 

Huxley et al 1996 
de Visser et al 1997 
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Table 4. A selection of miscellaneous studies 



Chemical formula 



Remarks 



Reference 



AMoeSs-xSez, A= Sn, Pb 
k-[BEDT-TTF] 2 Cu[NCS] 2 
Bi 

a-0 2 
Sb 

PrNi 5 
RbC 60 

Be 



Penetration depth measurements 

Penetration depth measurement and fluxon dynamics 

theoretical investigation of line broadening due to 
inhomogeneity of the Landau orbital magnetization 

characterization of the magnetic phase transition 
below the a-fS transition temperature 

investigation of the anomalous Knight shift 

observation of muon induced crystal field effects 

magnetic phase diagram investigation 



diamagnetic domains; see section 3.3 



Birrer et al 1993 
Harshman et al 1994 
Solt 1994 

Storchak et al 1994 

Birrer et al 1995 

Feyerherm et al 1995b 

Cristofolini et al 1995 
MacFarlane et al 1995 
Uemura et al 1995 

Solt et al 1996a 
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Table 5. A selection of work carried out on superconducting and non- 
superconducting oxides 



Chemical formula 


Remarks 




Reference 


Superconducting oxides 



Bi2+iSr 2 _ :!; CaCu20 8+5 

HgBa2 Ca3 Cu4 10+ ^ 

B12 Sr2 Cai _ x Y x Cu2 O g+ s 
Bi2_ a: Pb :l ;Sr2CaCu20 8+5 

YBaaCuaOy 

YBaaCusOs.e 

YBa2(Cu .96Zno.o4) 3 :E 

La2Cu04 
YBa 2 Cu 3 06.95 

La2-iSrj;Cu04 

La 2 Cui_ a; Zn a ;04 

YBa 2 Cu 4 8 H :E 
YBa^CuyOz-a 

YBa 2 Cu 3 H 

Yo. 8 Cao.2Ba2(Cui_ y Zn y ) 3 07_ 

La 2 _ a: Sr a; Ni04 + i 
Nd 2 _ :c Ce ;c Cu04 

YBaaCCui.yZnyJgO,; 

La2- a: Sr I (Cui_j / Zn H )04 

La2Cu04+i/ 



flux-line lattice study : po ssible observation 
of melting; see section p. 2| 

test of the univ ersal behaviour; 
see section |5.l| 

results unconsistent with universal behaviour; 



see section 5.1 



theoretical analysis of the damping of the 
precession component 

investigation of effect of sulfur substitution 
through measurement of the penetration depth 

mapping of the magnetic phase diagram for 
6.00 < x < 6.92 

computation of the hyperfine field 

study of the flux-line lattice field distribu tion 
for a mosaic of single crystals; see section |5.l| 

NMR and fiSK study of the weakly doped region; 
possible evidence for phase separation 

NMR and fiSK study of the effects of substitution 
of magnetic Cu 2+ with diamagnetic Zn 2+ 

study of the antiferromagnetic order 

in-plane London penetration depth anisotr opic due 
to superconductivity in chains; see section p,l| 

sintered samples; unconsistent with results of Sonier 
et al 1994 and 1997 and Riseman et al 1995 
for y = 6.95; see section p.l| 

dependence of the condensate density on Zmdoping 
points towards d-wscve pairing; see section p.l| 

magnetic phase diagram investigation 

no magnetic order but spin fluctuations for x = 0.02; 
ordering of the Nd moments at low T for x = 

effect of Zn doping on the superconducting electron 
density; see section p,l| 

investigation of the spin-glass state 



Lee et al 1993, 1995 and 1997 
Bcrnhard et al 1995a 
Aegerter et al 1996 

Niedermayer et al 1993 
Weber et al 1993 
Aristov and Maleyev 1994 
Cloots et al 1994 
Mendels et al 1994 

Shukri B Sulaiman et al 1994 

Sonier et al 1994 and 1997 
Riseman et al 1995 

Borsa et al 1995 
Corti et al 1995 



Gliickler et al 1995 

Tallon et al 1995 and 1996 
Bernhard et al 1995b 

Zimmermann et al 1995 



Bcrnhard et al 1996 

Chow et al 1996 
Hillberg et al 1997 

Nachumi et al 1996 

Pomjakushin et al 1996 
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Table 6. A selection of work carried out on superconducting and non- 
superconducting oxides (continuation of Table Bh 



Chemical formula 



Comments 



References 



Non superconducting oxides 



Cao.86Sro.i4Cu02, Sr2Cu03 
CuGe03 and Cui-^Zn^GeOs 



SrCr 8 Ga 4 Oi 9 

Y 2 Cu 2 5 
RNi0 3 

Sr„_iCu n+ i02n 

(Y 2 - !C Ca !C )Ba(Nii_j / Mg !/ )0 5 

La 2 Co a; Cui_ ;c 04 +a 
Y2M02O7, Tb22Mo 2 7 

CuO 

LaCu02.5 

KTr 3 (OH) 6 (S0 4 ) 2 
Tr = Cr, Fe 

La2Cui_ a; Li :E 04 

Ca 2 Cu0 3 , Sr 2 Cu0 3 

LiV 2 4 

SrCu0 2 



magnetic order in infinite-layer and chain 
possible spin-glass ground state of the doped 
compound; see comments of Kadono 1997 



a spin-liquid ground state is proposed; 
see comments of Dunsiger et al 1996 

magnetic phase diagram investigation 

characterization of the magnetic order in the 
low-temperature insulating phase 

magnetic phase diag ram of spin ladder 
systems; see section |j.2| 

Haldane system; chain length controlled 

by doping; effect of doping on the ground state 

magnetic phase diagram investigation 

spin dynamics of geometrically frustrated 
magnets; large density of states near zero energy 

magnetic phase diagram investigation 

the ground state is magnetically ordered 
rather than spin-liquid 

long range ordering in the Fe compound; 
no such ordering in the Cr compound 

magnetic phase diagram; formation of a singlet 
ground state at large doping 



infinite-chain cuprates; see section 3.2 

heavy fermion oxide with no magnetic ordering 

magnetic phase transition at ~ 2 K 
in this zigzag chain compound 



Keren et al 1993 and 1995 
Lappas et al 1994 
Garcia-Mufioz et al 1995b 
Sohma et al 1995 
Tchernyshyov et al 1995 

Uemura et al 1994 

Duginov et al 1995 
Garcia-Mufioz et al 1995a 

Kojima et al 1995a 

Kojima et al 1995b 

Lappas et al 1995 
Dunsiger et al 1996 

Grebinnik et al 1996 
Kadono et al 1996 

Keren et al 1996a 

Le et al 1996b 

Kojima et al 1997 
Kondo et al 1997 
Matsuda et al 1997 
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Appendix A. The magnetic field at the muon site 

We first express the measured local magnetic field in terms of the field components in 
the crystal reference frame. Then we write these field components in terms of their 
spatial Fourier components. Finally, we complete this appendix by providing a method 
of computing the field distribution at the muon site for a given magnetic structure. 

Appendix A.l. The magnetic field in the laboratory and crystal reference frames 

Whereas we are interested in the characterization of the magnetic field components at 
the muon site in terms of parameters of the crystal under study, the measurements are 
done in the laboratory reference frame. Here we express the linear relation between 
the field components in the two frames. 

For simplicity we suppose that the magnet has only one type of localized magnetic 
moment and that the muon occupies only one interstitial site. We must consider two 
orthonormal reference frames: the laboratory reference frame (X, Y, Z) where X, Y, 
and Z are unit vectors and a reference frame (x, y, z) attached to the crystal axes. Its 
unit vectors are chosen parallel to the crystal axes according to the symmetry of the 
compound. (X, Y, Z) are defined in the (x, y, z) frame through the Euler angles 9, p 
and tp. Note that 9 and tp are also the polar angles of the Z axis in the (x, y, z) frame. 
The components of the local field in the laboratory frame are written as a function of 
the components in the crystal frame: Bf oc = ^ a R ea (0, tp, VOAoc- We ^ ne n ' ne 
R ga components: 

Rxx = cos tp cos tp cos 9 — sin tp sin tp, R xy = cos tp sin tp cos 9 + sin tp cos tp, ( Al ) 



Ryx = — sin tp cos tp cos 9 — cos tp sin tp, R yy = — sin tp sin p cos 9 + cos tp cos tp, (A2) 



Rxz = — cos tp sin 9, R yz = sin tp sin 9, (A3) 



Rzx — cos tp sin 9, R zy = sin p sin 9, R zz — cos 9. (A4) 

Since R is a unitary matrix, it is easy to express f?g c in terms of Bf oc using the 
matrix elements listed above. 

Appendix A. 2. The Fourier components of the magnetic field 

This appendix is based on an extension of the work of Yaouanc et al 1993b. The 
reader will consult with profit the paper of Lovesey and Engdahl 1995. 

We express the local field in the crystal frame as a function, on one hand, of a 
tensor which describes the coupling between the localized spins of the magnet and the 
muon spin and, on the other hand, of the localized spins components themselves. The 
a component of the field is given by 




(A5) 
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v is the volume per magnetic atom, rid the number of magnetic atoms in the cell used 
for the description of the magnet, ql the Lande factor of spin Ji^d which is at the 
distance vector Yi + d from the muon, \ib the Bohr magneton, (j,q the permeability of 
free space and {a, (3}={x, y, z}. r defines the location of the muon relative to the 



origin of the crystal lattice; see figure Al. Note that the equivalent crystallographic 



sites allowed by the point group symmetry at the muon may give rise to different 
G^(q) and therefore to inequivalent magnetic sites, each characterized by a given r . 
The index i runs over the cells and d over the sites inside a cell. Since one of the most 
important characteristics of a magnetic structure is its periodicity we introduce the 



spatial- Fourier component J^(q). Using the vectors defined in figure Al we derive 



£ ££Gf(q)e X pHq-r )jf(q), (A6) 

*VI lid" 

p—x,y : z d q 

where we have defined 

Gf (q)=£G?£„exp(iq.r i+( ,) 1 (A7) 



and 



^(q) = £-C ex pHq-(i + d)]. (as) 




Figure Al. Definition of vectors relative to the crystallographic (+, x) and muon 
(□) sites. The points (•) specify the origin of each cell. The point at the extremity 
of ro is the origin of the crystal lattice. The drawing is done for a two dimensional 
square lattice with two magnetic atoms per cell. See the text for the definitions of 
the different vectors. 

The tensor G^ /3 (q) is the sum of two tensors: G^(q) = D^(q) + H^(q). 
Whereas D^(q) describes the dipolar interaction between the spins of the magnet 
and the muon spin, H^(q) accounts for their hyperfine interaction. It is convenient 
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to rewrite 



An Pf (q) 



cf (q) 



(A9) 



where P^^(q) is the longitudinal projection operator and C^(q) a symmetric tensor 
which reflects the point group symmetry at the muon site. The expression of C^(q) 
is computed using Ewald's method: 



4p 
K K^O 

n d vg 3 \ - 

^(tt)*/*^ 



1 — exp 



4g 2 

f 9/3) ¥>o 



(q + K) ; 



2 r l+d , Q r. J+di/ 3^3 /2 (<?V 2 +d ) - <5 Q/3 - ' ■-•• 2 



exp (-iK • r 0+t z) 

{Q' r i+d)\ ex P (*q • r «+<i) » ( A1 °) 



where we use the Misra functions: 

/oo 
df3(3 m exp(-/9x). 



(All) 



K is a vector of the reciprocal lattice. Expression (A10) gives the same result for all 
values of the Ewald parameter p, but for numerical applications a value of g is chosen 
which ensures that both series of (A10) converge rapidly. 

Whereas P^(q) is only piecewise continous at q = 0, reflecting the long range 
nature of the dipolar interaction, C?^(q) is analytical. The presence of the projection 
operator in (A9) is important for understanding the critical behaviour of the relaxation 
rate in ferromagne ts a nd the temperature dependence of this rate in the spin wave 
regime (see se ctio n 4.1). 

In section 2.3 we have noted that Bi oc can be written as a sum of seven terms. We 
can identify the first five terms. (B con + Bt ra ns) results from the polarized conduction 
electrons and the transferred fields. To first approximation, these two fields are 
isotropic in metals. The (B con + Btrans) contribution to the local field is obtained by 
substituting H for G in (A6). In the same way, the (B^j +Bl) and Bd em contributions 
are derived from (A6) by substituting 47rC and —AttPl for G, respectively. 

Although the spatial-Fourier transform may seem a complicated procedure for a 
result which can be obtained by a simple lattice sum, Appendix A. 3 and Appendix C 
will show its unique possibilities. 



Appendix A. 3. The field distribution at the muon site 

Even for a simple magnetic structure, crystallographically equivalent sites may not be 
magnetically equivalent, i.e. more than one local magnetic field may exist. This is 
nicely explained for the case of Fe by Seeger 1978. If the number of fields is sufficiently 
large and their values are close, the muons will probe a field distribution which may be 
far from Gaussian as supposed in section || The relation between a magnetic model 
and the field distribution at the muon site, £>(Bi oc ), was first investigated by Szeto 
1987. 

In an experiment the several million implanted muons may enter different crystal 
sites which are related by the lattice translation periodicity. Therefore an average 
over these r sites is required: r = ro + Zx + my + nz where {Z, m, n} are integers. We 
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recall here that if the point group symmetry of the muon site gives several equivalent 
positions, one has to consider each of these as they may give rise to different coupling 
tensors G<j(q) i.e. to non equivalent magnetic sites. 

For the sake of simplicity, we consider a magnetic structure described by a single 
propagation vector k. This represents the majority of structures found in nature, and 
the extension to multi-k structures is easily performed. Then the Fourier components 
of the magnetic moment may be written as 

M d (q) = M d (pk)<5(q - pk), (A12) 

p 

where 5 is the Kronecker symbol. The sum is over the harmonics of the magnetic 
structure; p = ±1 for a sinusoidal modulation. Note that M<j(q) is related to Jd(q) 
via the relation M<;(q) = gLfJ,B^d(<l)- Using (A12) and (A6), we derive 

^ r ) = £^££ E Gf(pk)M^k)ex P Hpk.r) (A13) 

d p f3—x,y,z 

For a commensurate magnetic structure k = (r/s)G where r and s are integers with 
r < s and G a reciprocal lattice vector. From the well known properties of G, we 
derive k-r = k-ro + 2ir(r/s)u where u is an integer. Using this result, the exponential 
term of (A13) is rewritten as 



exp[— ipk ■ r] = exp(— ipk ■ r ) 



2irpr 
cxp I i 



(A14) 



During a measurement, an ensemble average is made, so that the integer u takes an 
enormous number of values. But, since {exp [i(2npr)/s)]} u as a function of u is a 
periodic function of period s, this phase factor takes a maximum of s different values. 

Since B^ c (r) depends linearly on B^ oc (r) and this latter field is a sum of Fourier 
components, Bi oc (r) is a sum Fourier components. Taking into account the property 
of the phase factor, we deduce that each component can take many different values. 

We discuss this important result. For simplicity we consider a primitive Bravais 
lattice (rid = 1), a sinusoidal modulation and suppose that the laboratory and crystal 
frames are aligned with Z and z parallel. For instance, if r = 1 and s = 3, up to three 
different fields can exist. If s is big, a large number of fields can result and in the 
limit, a quasi-continuous distribution is generated. 

To proceed further analytically wc consider this limit and write Bi oc (r) = 
cos[27rw(r)]B m with < w < 1. This form is strictly valid for an incommensurate 
magnetic structure. Using the well known formula for a distribution 

D (Bioc) = ((S [Bioc - Bi oc (r)])) , (A15) 

where ((...)) stands for the spatial average over r, we derive (Le et al 1993a) 

D(Bta, -WnE' (A16) 

for < Bioc < B m and D (Bi oc ) = otherwise. If we suppose that Bi oc is 
perpendicular to Z, we find using (3) 

Pz(t) = Jo(^B m t), (A17) 
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where J is a Bessel function. When ^^B m t is large, the model predicts Pz{t) = 
^/2/ (777^ i? TO t) cos(7 jU B TO t — 7r/4) instead of the usual Pz(t) = cos(7 Al i? m t), i.e. the 
depolarization function presents at large t a damped oscillation shifted by 45°. 

If the muon site is known, using the B^ c (r) expression (A13) it is possible to 
compute numerically the field distribution for any magnetic structure. While with the 
/iSR technique it is not possible to determine in detail a magnetic structure, one can 
test proposed structures. 



Appendix B. Muon spin relaxation in a longitudinal field 

When the dynamics of the magnetic field at the muon site is sufficiently rapid, Pz(t) 
takes the form (McMullen and Zaremba 1978 and Dalmas de Reotier and Yaouanc 
1992) 

Pz(t)=exp[-ip z (t)], (B18) 



with 



i>z{t)=ll J dr (t-r){cosM[$ M (r) + ^(T)] 

+ sin(^r) [$ xy (r) - <P yx (t)] } . (B19) 
3>(t) is the symmetrized correlation-tensor of the magnetic field at the muon site: 

$°*V) = \ {(Kc(t)bQ + (< c K c (r))] (B20) 

(...) stands for the thermal average. 

If the fluctuations arc sufficiently rapid, one neglects r in the (t — r) factor and 
extends the integral to infinity. Then Pz(t) is an exponential function characterized 
by the relaxation rate Xz — ipz{t)/t. 

In zero-field $ Q/3 (r) is an even function of t. In an applied field, this may not be 
true because time reversal symmetry is broken. However, in many cases the breaking 
terms in the Hamiltonian are so small that the effect of the field is negligible. Then, 
even in an applied field, $ a/3 (r) is an even function of t. We will suppose that this 
property holds. This enables the integration over r in (B19) to be extended from —00 
to 00. 

In zero field Xz is given by 

Az = 4f/ dr[^ xx (T) + ^ YY (r)]. (B21) 

^ J-00 

This simple result can be derived from the Fermi golden rule (Lovesey et al 1992). 

In a longitudinal field, a term proportional to sin(w At r) is present as seen in (B19). 
In most cases this term is zero either because time reversal symmetry is not broken 
or uj^Tc <C 1 where t c is a characteristic time of the fluctuations. Then we have 

X z = ^ dr cos Kr) [$ xx (r) + <f> YY (t)] . (B22) 

Z J -co 

As stressed by Dalmas de Reotier et al 1996, Pz(t) is, in general, an exponential 
function only if no spatial average of the depolarization function is needed. This means 
that Pz(t) has, in general, no reason to be an exponential function for measurements 
on a polycrystalline sample. 



49 



Appendix C. Longitudinal relaxation rate and correlation-functions 

In this section we analyze the relaxation rate in terms of the spin-spin correlation- 
tensor of the magnet. 

Using the results presented in the two previous sections, $(r) can be written in 
terms of correlation- functions of the spatial Fourier components of the lattice spins. 
It is useful to introduce the time- Fourier transform of a function /(t): 

/M = h I dT ex p(-^ T )/( r )- ( C23 ) 

Using (A6) and (B22), Xz can be expressed in terms of the spin-spin correlation-tensor 
between spins belonging to sublattices d and d', Ad,d'(q, u>). With the definition 

A£,'(q,w) = \ [(jJ(q,c)Jj:(-q)) + ( J^(-q) Jj(q, , (C24) 
the Xz expression writes 

d 3 q 



= ^J2L 0a (9,<p) J ^EE^^N^fe^). (C25) 

(3, a 7,7' d 7 d' 



The sum over q in (A6) has been replaced by an integral over the first Brillouin zone 
following the usual rule J2q ~^ I ^/(2"") 3 d 3 q where V is the volume of the sample. 
We define V — (^ /47r) 2 7 2 (g^/Lis) 2 . L(9,ip) is a symmetric matrix which accounts 
for the fact that the symmetry axes of the magnet do not necessarily coincide with 



the laboratory frame axes. 9 and ip have been defined in Appendix A.l. We have 



L xx = cos 2 tp cos 2 9 + sin 2 ip, L yy = sin 2 ip cos 2 + cos 2 <p, L zz — sin 2 9, (C26) 

L X y = — cos ip sin ip sin 2 9, L xz = — cos <p cos sin 0, L yz = — sin </? cos 9 sin 0. (C27) 

It is often useful to consider A(q) = J^d d' A<j,d'(<j) an d the average coupling tensor 
G(q) = l/iiiErf^f !) since for some non Bravais crystal structures such as the 
hexagonal closed compact structure (Dalmas de Reotier and Yaouanc 1994) (C25) 
can be written in terms of A(q) and G(q), i.e. the sum over d and d! disappears in 
the expression of Xz- 

The symmetry properties at the muon site and of the magnet itself leads to 
considerable simplification. Examples are found in Yaouanc et al 1993a and 1993b, 
Dalmas de Reotier and Yaouanc 1994, Bonville et al 1996, Dalmas de Reotier et al 



1996, Yaouanc et al 1996a and in section 4.1 



As an example of the drastic simplifications which can occur, we suppose that only 
one spin-spin correlation matters. This is the case for an Ising system or a magnet 
with its magnetic ions in a cubic environment. Then the following formula can be 
derived: 

f d 3 a n f°° drcos(wur)(J 2 (q,r)J z (-q)) 
X z = v -^-A 2 (q)^ \ M A ™ ' 1 . (C28) 



(2tt) 3 ™ J(J + l)/3 

Xz depends on the function A 2 (q) which accounts for the coupling between the muon 
and the lattice spins. In zero-field, when the intersite correlations are neglected (A 2 (q) 
= A 2 <5(q)) and the spin-spin-correlation is an exponential function characterized by a 
fluctuation rate v, we recover the motional narrowing result: Xz = 2A 2 /v, see (10). 
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Appendix D. The /iSR response function for hard superconductors 

In this appendix we calculate the magnetic field inside the mixed state of a hard 
superconductor. For such a superconductor the Ginzburg-Landau parameter k is 
much larger than 1 (k = A/£ where A is the penetration length and £ the coherence 
length) . 

It is well known that a vortex lattice is formed in the mixed state. We suppose 
that this lattice is ideal. This means we disregard pinning and vortex "phases" such 
as the glassy or liquid states (Blatter et al 1994, Brandt 1995). Pinning can strongly 
influence the results of a /xSR experiment as shown in details by Wu et al 1993. We 
consider the case where the field B cxt is parallel to one of the three main axes a, b 
and c of the penetration-length tensor. With these simple conditions, the vortices 
are straight field tubes parallel to B cxt and form a regular lattice. We define an 
orthonormal reference frame (x,y,z) such that the vortex tubes are also along z. We 
denote A a , Ah and A c , the penetration lengths for currents flowing along the a, b and 
c axes, respectively. 

As shown by Kogan (Kogan 1981), for a conventional vortex lattice and B cxt 3> B cl 
where B c i is the lower critical field, the angle characterizing the lattice depends only 
on the penetration-length ratio : 

tana,, = V3(X x /X y ). (D29) 

As expected, if the penetration length is isotropic, a v — 60°. Therefore, although 
a v is most naturally measured by small angle neutron scattering (SANS), one can 



determine this angle by /iSR. An example is given in section 5.3 
It is convenient to introduce the Fourier components 

B(G) = J B(r) exp(—iG ■ r)d 2 r/ S (D30) 

of the periodic magnetic field 

B(r) = ^B(G)exp(iG-r), (D31) 
G 

where G are the vectors of the vortex reciprocal lattice and S the surface of the vortex 
lattice unit cell. Since B X (G) = B y (G) — 0, there is no transverse field component. 
Interestingly, the only non-zero component, B Z (G), is the form factor measured by 
SANS. 

Using Kogan's formula, the form factor factorizes (Yaouanc et al 1997a), B z (G pq ) 
= Bq ■ c pq (b), where 

1 CiV' 2 $0 



S » = ^UJ xS~, (D32) 

c pq {b) are universal functions of b = B / B& where B is the mean induction, which for 
26k 2 > 1 may be equated to B cx t- {p, q} are the two indices denoting a Bragg peak. 
Recently the c pq (b) functions have been computed numerically in the framework of 
the conventional Ginzburg-Landau theory (Brandt 1997). In terms of the functions 
b pq (b) defined by Brandt, we have c pq {b) = b pq (b)/(p 2 — pq + q 2 ). 

Since a typical penetration length is much larger than a crystal lattice parameter 
which is the characteristic distance between adjacent muon stopping sites, the muon is 
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a good probe for the vortex field distribution. Then it is straightforward to compute 
the distribution using (A15). This has been done by many authors (for example, 
Sonier et al 1994 and Greer and Kossler 1995). A characteristic of the flux line 
lattice observed in many experiments (for example, see figure [l8] and figure ^0|) is 
the pronounced tail towards high fields due to the vortex cores. A useful method 
to characterize a distribution is to consider its moments or, since (B z ) ^ 0, A™ = 
((B z — (B z )) n ). The variance which is a measure of the width of the distribution 
separates into two factors, -y/A 2 = Aq • f v {b) where 



A = 0.06092 



A..,- A,, 



(D33) 



is the London limit (£ x , £ y — > 0) (Barford and Gunn 1988) and f v (b) a universal 
function which accounts for the core size. <3?n = 2.07 x 10 -15 Tm 2 is the quantum of 
flux and the prefactor 0.06092 is a pure number. f v {b) computed using the data of 
Brandt 1997 is plotted in figure A2. It is strongly field dependent even at low reduced 
fields b. The shape of a distribution is characterized by its skewness parameter : a 
= (Af) 1 / 3 / (A 2 ,) 1 / 2 . For a symmetric distribution a = 0. a(b) is also presented in 
figure |A2| . Note the strong field dependence of a near 6 = 0. We have a(b = 0) = 
1.446. 



Dalmas de Reotier and Yaouanc; fig. A2a 



Dalmas de Reotier and Yaouanc; fig. A2b 




0.0 0.2 




Figure A2. Field dependence of the functions f v and a. b = B/B c 2 where B is 
the mean induction, which for 2bK 2 > 1 may be equated to -B C xt- In the inserts we 
present the functions versus fr 1 / 2 to extend the low b part, (adapted from Yaouanc 
et al 1997a and Brandt et al 1997). 



The results presented in figure A2 have been deduced from the numerical solution 
of the conventional Ginzburg-Landau expansion (Brandt 1997). Although this 
expansion has its own limit (basically it is valid close to the superconducting phase 
transition in the temperature-field phase diagram), a comparison between the results 
derived from the Ginzburg-Landau model and the commonly used Gaussian cutoff 
London model shows that the latter model strongly underestimates the effect of the 
field on the variance (Yaouanc et al 1997a). In fact, the Gaussian model does not 
have any theoretical support. 

The motion of the vortices on a time scale shorter than a few microsecond may 
significantly effect the field distribution as probed by muons. Since the amplitude of 
such fluctuations is field dependent, the variance can exhibits a strong field dependence 
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as shown by Song 1995. 

The Ginzburg-Landau model does not include the possibility of a Knight shift. 
Because this shift is expected to be different in the superconducting and in the normal 
regions, and since normal (vortex core) and superconducting regions coexist in the 
mixed state, the measured field distribution is a convolution of the distributions due 
to the Knight shift and the flux line lattice. This point was first raised by Feyerherm 
et al 1994. 

Since a SANS experiment measures |S 2 (G)|, one may estimate that SANS and 
/iSR experiments measure the same physical quantity. This idea seems to be supported 
by the fact that c w and f v have approximately the same field dependence (Yaouanc et 
al 1997a). But there is at least one important difference. In order to observe a Bragg 
peak by SANS, the correlation length of the vortex lattice must be sufficient. This is 
not required to observe a vortex field distribution by fiSR which is a local technique. 

A "pure" field distribution is not observed because of broadening effects due to 
the finite lifetime of the muon, the defects of the vortex lattice, the demagnetization 
field and the nuclear dipolar broadening. We refer the reader to Riseman et al 1995 
and Schneider et al 1995 for more information. 

Up to now, we have only considered a three dimensional vortex lattice. But when 
the anisotropy of the penetration length is strong, it is more appropriate to describe 
the vortices as two dimensional "pancakes" vortices in individual coupled layers. 
Experimentally, Bi 2 + I Sr 2 - 2: CaCu208+ ( 5 corresponds to this case (see section ^2). As 
first noticed by Brandt 1991, the misalignment between the layers leads to a variance 
smaller than for a standard three dimensional flux line lattice. A comparison between 
simulated field distributions for three and two dimensional vortex lattices has been 
made by Schneider et al 1995. 

One of the most interesting parameters deduced from the investigation of a 3D 
field distribution is the penetration length A; indeed its determination is one way to 
probe the nature of the low energy excitations and therefore the pairing state of a 
superconductor. The relation between A and these excitations is reviewed by Gross et 
al 1986 and Gross- Alltag et al 1991. 
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